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PROCESSING WAVEFORM-BASED NDE 

This application is a continuation-in-part of application 
Ser. No. 12/163,382 filed Jun. 27, 2008 now U.S. Pat. No. 
7,876,423. The instant application and application Ser. No. 5 
12/163 ,3 82 filed Jun. 27, 2008 have the same inventor and are 
commonly owned by the same assignee. 

ORIGIN OF THE INVENTION 

to 

The invention described herein was made by an employee 
of the United States Government, and may be manufactured 
and used by the government for government purposes without 
the payment of any royalties therein and therefor. 

15 

FIELD OF THE INVENTION 

As a result of Space Shuttle Columbia Accident Investiga- 
tion Board recommendations, an aggressive program to 
eliminate all External Tank Thermal Protection System 20 
debris-shedding at the source was initiated. Terahertz c-scan 
imaging is an emerging and very effective nondestructive 
evaluation (NDE) technique used for dielectric materials 
analysis and quality control in the pharmaceutical, biomedi- 
cal, security, materials characterization, and aerospace indus- 25 
tries. 

BACKGROUND OF THE INVENTION 

Flaws present in the Space Shuttle external tank thermal 30 
protection system may play a role in foam release and are 
therefore important to detect and characterize prior to flight. 
The external tank configuration has sprayed-on foam insula- 
tion placed on top of the metal container and thus lends itself 
to terahertz inspection. Terahertz inspection has shown sig- 35 
nificant promise for detection of voids in the foam. Other 
potentially undesirable foam anomalies that have been iden- 
tified by NASA include density variations and crushed foam. 
Velocity imaging can be used to identify density variations. 

Terahertz waves are electromagnetic waves with wave- 40 
lengths on the order of 200 to 1000 pm. Reflections occur to 
varying degrees at interfaces between materials with dissimi- 
lar dielectric properties (difference in indices of refraction). 
Metallic materials totally reflect terahertz waves while non- 
polar liquids, dielectric solids, and gases are at least partially 45 
transparent to terahertz energy. Continuous wave (narrow- 
band) and pulsed (broadband) terahertz systems exist. 

Several attempts to separate thickness and micro structural 
variation effects in ultrasonic images are noted in the litera- 
ture. Several references showed single point (non-imaging) 50 
ultrasonic measurement methodology that accounted for 
thickness variation effects. See, Sollish. B. D., Ultrasonic 
Velocity and Thickness Gage, U.S. Pat. No. 4,056,970, Nov. 

8. 1977. Hsu, D. K. et al., Simultaneous determination of 
ultrasonic velocity, plate thickness and wedge angle using 55 
one-sided contact measurements, NDT&E International 
1994 vol. 27, no. 2, pp. 75-82 and Piche, L., Ultrasonic 
velocity measurement for the determination of density in 
polyethylene, Polymer Engineering and Science, vol. 24, no. 

17, Mid-December 1984 pp. 1354-1358. Hsu et. al, 1994, 60 
simultaneously determined ultrasonic velocity, plate thick- 
ness and wedge angle. Piche, 1984, described a single point 
ultrasonic velocity measurement method using a reflector 
plate located behind the sample that does not require prior 
knowledge of sample thickness and lends itself to multiple 65 
measurements within a sample of nonuniform thickness. Sev- 
eral references proceeded to scale up and automate this ultra- 
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sonic method to obtain ultrasonic velocity images for plate 
and cylindrical samples of various materials of non-uniform 
thickness. See, for example, Dayal, V., ‘An Automated 
Simultaneous Measurement of Thickness and Wave Speed by 
Ultrasound,” Experimental Mechanics, 32(3), pp. 197-202, 
1992; and, Roth, D. J., Carney, D. V., Baaklini, G. Y., Bodis, 
James R., Rauser, Richard W., “A Novel Ultrasonic Method 
for Characterizing Micro structural Gradients in Tubular 
Structures,” Materials Evaluation, Vol. 56, No. 9, September 
1998, pp. 1053-1061. 

A procedure utilized in ultrasonics and terahertz in which 
the substrate reflector plate time-of-flight scan with no 
sample present is subtracted from the same scan with the 
sample in place is useful to characterize microstructure and 
correct for setup nonuniformity i.e., levelness, but it will not 
separate thickness and micro structural effects. 

Ultrasonic methods to simultaneously measure or charac- 
terize thickness and density (or variation as such) require 
water coupling. Additionally, the ultrasonic methods cannot 
be used for foam inspections due to the highly porous nature 
or highly cellular structure of foams. The terahertz method is 
totally non-contact, requires no coupling, and works in air. 

SUMMARY OF THE INVENTION 

Terahertz imaging is being used at NASA for nondestruc- 
tive evaluation of the Space Shuttle external tank thermal 
protection system sprayed-on foam insulation (SOFI). The 
NASA Engineering and Safety Center tasked a technical team 
to develop improved inspection methods to characterize foam 
anomalies to help alleviate foam shedding on the space 
shuttle tanks. Foam density variation was identified as a 
potential problem in which thermal expansion mismatch 
between areas of different density could result in crack for- 
mation, subsequent foam shedding, and endangerment of the 
space shuttle orbiter. Prior to implementation of the method 
disclosed herein it was not possible to quantitatively measure 
density using a totally non-contact, non-water-coupled 
method. Generally, terahertz is used in the pulse-echo c-scan 
configuration to map variations in the peak amplitude of the 
echo off of the metal substrate (equivalent to the location of 
the back surface of the foam) that occur when scanning across 
a section of foam in order to detect voids, cracks, disbonds, 
and any sort of discontinuity. Traditional c-scan imaging 
scales the peak amplitude values (to an 8- or 16-bit gray or 
color scale) at each scan location to form an image. 

Since the pulse-echo terahertz method results in a wave- 
form with echos being received off of the front surface of a 
dielectric material and a metal (electrically-conducting) sub- 
strate that the dielectric material rests on, obtaining the time 
delay between front surface and substrate (with the sample 
present) echos is possible. Terahertz velocity is affected by 
variations in a volumetric micro structural property such as 
physical density and thereby once a relationship between the 
two variables is established, a non-contact precise measure- 
ment of density can be made using terahertz energy. 

If the dielectric material has flat and parallel sides such that 
no thickness variation exists, the time delay between the front 
surface echo and substrate echo with the sample present will 
be indicative of only micro structural variation. By obtaining 
the relationship between velocity and a microstructure prop- 
erty, such as density, using a series of samples of different 
density, one can then predict the density of the material and 
subsequently map density variations within the material 
using the established relation between velocity and density. 

A pulse-echo terahertz velocity measurement is made by 
sending terahertz energy via a transceiver (device that has 
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both a transmitter and a receiver) into and through a dielectric 
(insulating) material (such as the shuttle external tank thermal 
protection system sprayed on insulating foam) backed by a 
metallic (electrically -conducting) plate that reflects the tera- 
hertz energy back to the transceiver. The terahertz transceiver 
is separated from the dielectric sample by an air path. Velocity 
(V) values are calculated using the time delay between the 
front surface echo (FS) and substrate/reflector plate echo 
(BS). With a dielectric sample present between the trans- 
ceiver and the reflector plate, the pulse that travels from the 
transceiver through the sample to the reflector plate (equiva- 
lent to the sample back surface position) and back to the 
transceiver is labeled BS and will be observed at time t\ Thus 
two “echos,” FS and BS, can have their peak positions in time 
measured and the time difference or time delay between them 
is determined. Alternatively, the entire echos may be cross- 
correlated to obtain the precise time delay between them. If 
thickness is non-variable in the sample, the time-of-flight 
and/or velocity measured will be indicative only of the micro- 
structure. The FS echo may require specialized signal pro- 
cessing to denoise and amplify it. 

A process for measuring the velocity of terahertz electro- 
magnetic radiation in a material sample without prior knowl- 
edge of the thickness of the sample is disclosed and claimed. 
Terahertz electromagnetic radiation is produced by a source 
(transceiver) spaced apart from the substrate and propagated 
at the speed of light, c, in a medium located between the 
source (transceiver) and the substrate. The process for mea- 
suring the velocity (independently of thickness) includes the 
following steps: emitting terahertz electromagnetic radiation 
from the source; measuring the travel time of the terahertz 
electromagnetic radiation to and from the substrate without 
the sample present, t" (t-double prime); placing the sample on 
the substrate; measuring the travel time of the terahertz elec- 
tromagnetic radiation to and from the substrate with the 
sample present, f (t -prime); subtracting t" fromt' to determine 
the transmission time difference with the sample present and 
without the sample present, At; measuring the travel time of 
the terahertz electromagnetic radiation to and from the 
sample, 2 time determining the difference between the FS 
echo and the BS echo, 2 t; dividing At by 2x and determining 
the quotient; and, subtracting the quotient from 1 to obtain a 
factor; multiplying the factor by, c, to determine the velocity 
of the terahertz electromagnetic radiation in the sample. The 
further step of determining the micro structural variation of 
the sample according to an algorithm is performed to deter- 
mine, for instance, the density of the foam. The step of deter- 
mining the micro structual variation of the sample includes 
determining the density of the material. The sample may be 
foam or another dielectric such as silicon nitride. The radia- 
tion may be pulsed or it may be continuous. 

If the surface to be evaluated is large, then mapping micro- 
structural variations in a plurality of locations is performed. 
When the terahertz electromagnetic radiation reaches the 
front surface of the sample, the echo therefrom may not be 
very prominent if the dielectric mismatch between the air path 
and the sample and the dielectric itself is not substantial 
enough. The step of measuring the travel time of the terahertz 
electromagnetic radiation to and from the sample, 2t l option- 
ally includes placing a sheet of plastic transparent paper onto 
the sample to create an adequate dielectric mismatch between 
the air and the sample. The plastic transparent paper is thin, 
for example, it may be approximately 250 pm thick or less. 
The plastic paper does not have to be transparent. It is neces- 
sary to know, a priori, the approximate distance between the 
transceiver and the front surface of the sample. As used herein 
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“a priori” means before knowledge of the exact distances 
between the transceiver and the front surface of the sample. 

The front surface echo (FS) from the dielectric material 
(sample) may be of very low signal-to-noise ratio (SNR) 
5 depending on the dielectric match between air and the 
sample. If a good dielectric match exists, much of the tera- 
hertz energy will be transmitted into the sample. Additionally, 
the focal plane sensitivity of the terahertz method disclosed 
herein, may, for samples of nonuniform thickness, result in 
to the front surface echo (FS) too far out of focus and thus 
reduces the signal to noise ratio (SNR) even further, thus 
limiting the thickness variation over which the method can be 
used. The approximate time location of FS must be known “a 
priori” and the wavetrain examined manually through obser- 
15 vation of the signal on an oscilloscope trace to determine what 
special post-processing needs to be applied. In this way 
amplification and denoising the front surface signal (FS) can 
be achieved. Therefore, the further process steps of control- 
ling, approximately, the spacing between the source (trans- 
20 ceiver) and the substrate and the spacing between the source 
and the front surface are usually performed preliminarily if 
necessary. It is also necessary to know the approximate dis- 
tance between the transceiver and the metal substrate. 

The step of controlling the spacing between the source and 
25 the substrate includes determining, initially, the approximate 
time location of the sample from the terahertz radiation 
source as well as the approximate time location of the sub- 
strate from the terahertz radiation source. 

Identical scan data may be used for measuring the thick- 
30 ness of a material sample using terahertz electromagnetic 
radiation in a material sample without prior knowledge of the 
velocity of the terahertz electromagnetic radiation in the 
sample. In other words density and thickness can be obtained 
from the same scan data. The terahertz electromagnetic radia- 
35 tion is produced by a source (transceiver) spaced apart from 
the substrate and propagated at the speed of light, c, in a 
medium (usually air) located between the source and the 
substrate. The steps in the process include: emitting pulsed 
(or continuous) terahertz electromagnetic radiation from the 
40 source; measuring the travel time of the terahertz electromag- 
netic radiation to and from the substrate without the sample 
present, t" (t-double prime); placing the sample on the sub- 
strate; measuring the travel time of the terahertz electromag- 
netic radiation to and from the substrate with the sample 
45 present, t' (t-prime); subtracting t" from t' to determine the 
transmission time difference with the sample present and 
without the sample present, At; measuring the travel time of 
the terahertz electromagnetic radiation to and from the 
sample, 2t 1 ; determining the time difference between the FS 
50 echo and the BS echo, 2 t; subtracting At from 2 t to obtain a 
subtraction result; dividing the subtraction result by 2 to 
obtain a quotient; and, multiplying the quotient by, c, to 
obtain the thickness of the sample. Additionally, the method 
for determining thickness may include evaluating, in a plu- 
55 rality of locations, the sample for thickness variations and 
mapping the thickness variations by location. 

The inventor discloses and claims herein a process for 
simultaneously measuring the velocity of terahertz electro- 
magnetic radiation in a material sample without prior knowl- 
60 edge of the thickness of the sample and for measuring the 
thickness of a material sample using terahertz electromag- 
netic radiation in a material sample without prior knowledge 
of the velocity of the terahertz electromagnetic radiation in 
the sample. The pulsed terahertz electromagnetic radiation is 
65 produced by a source spaced apart from the sample under test 
and propagated at the speed of light, c, in a medium located 
between the source and the sample. The medium is typically 
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air. The process includes the steps of: emitting (pulsed or 
continuous) terahertz electromagnetic radiation from the 
source; measuring the travel time of the terahertz electromag- 
netic radiation to and from the substrate without the sample 
present, t" (t-double prime); placing the sample on the sub- 
strate; measuring the travel time of the terahertz electromag- 
netic radiation to and from the substrate with the sample 
present, f (t-prime); subtracting t" from t’ to determine the 
transmission time difference with the sample present and 
without the sample present, At; measuring the travel time of 
the terahertz electromagnetic radiation to and from the 
sample, 2 t l ; determining the time difference between the FS 
echo and the BS echo, 2 t; dividing At by 2x and determining 
the quotient; subtracting the quotient from 1 to obtain a factor; 
multiplying the factor by, c, to determine the velocity of the 
terahertz electromagnetic radiation in the sample; subtracting 
At from 2 t to obtain a subtraction result; dividing the subtrac- 
tion result by 2 to obtain a quotient; multiplying the quotient 
by, c, to obtain the thickness of the sample; and, evaluating, in 
a plurality of locations, the sample for micro structural varia- 
tions and for thickness variations, and mapping the micro- 
structural and thickness variations by location. 

Another process for simultaneously measuring the velocity 
of terahertz electromagnetic radiation in a material sample 
without prior knowledge of the thickness of the sample and 
for measuring the thickness of a material sample using tera- 
hertz electromagnetic radiation in a material sample without 
prior knowledge of the velocity of the terahertz electromag- 
netic radiation in the sample, the terahertz electromagnetic 
radiation produced by a source spaced apart from the sample 
under test and propagated at the speed of light, c, in a medium 
located between the source and the sample is disclosed herein 
which comprises the steps of: merging the FS, BS, and M" 
scan data sets into a fused data file; producing a “fused” 
waveform; calculating time delays between FS, BS and 
M M echos; determining precision thickness -independent 
velocity images that map micro structure; and, determining 
precision micro structure-independent thickness images that 
map thickness. The steps of calculating the time delays and 
determining precision thickness and micro structure are per- 
formed by cross-correlating FS and BS signals (echos) and by 
cross correlating BS and M" signals (echos). The steps of 
calculating the time delays and determining precision thick- 
ness and micro structure may also be performed by precisely 
identifying peaks of FS, BS and M". Where FS is believed to 
vary across a sample, an additional step of processing and 
gating the FS signal prior to merging the data into a fused data 
file is performed. Optionally, the step of amplifying the pro- 
cessed FS signal prior to merging the data into a fused data file 
is performed. To obtain a good dielectric mismatch between 
the sample and the air, the step of applying a dielectric sheet 
on the FS, followed by gating and processing the FS signal 
prior to merging the data into a fused data file is performed. 

The methodology disclosed herein has applicability to all 
dielectric materials where non-contact, non-water-immer- 
sion precision determination of microstructural (density) 
variation is required. It can be used for precision density 
mapping in dielectric ceramic materials, other types of foam, 
and dielectric composite materials. 

C-scan imaging involves mapping variations in the time- 
of-flight of a terahertz echo peak, or mapping the time delay 
between front surface and substrate (with the sample present) 
echos (FS, BS). The novel implementation described herein 
concerns itself more with mapping thickness or global micro- 
structural variation (such as physical density variation) as 
opposed to discrete flaw detection. Time delay between the 
front surface echo (FS) and substrate echo (BS) (with the 
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sample present) is directly affected by thickness variation (d) 
and terahertz velocity in the material (V). Terahertz velocity 
is affected by variations in a volumetric microstructural prop- 
erty such as physical density. 

5 The terahertz method of inspecting metal reflector-backed 
dielectric materials provides velocity images free of thick- 
ness variation effects, i.e. thickness -independent. Addition- 
ally, the same methodology can be slightly manipulated to 
obtain thickness images free of microstructural variation 
10 effects, i.e. microstructure-independent. In simple terms 
thickness can be measured without knowing velocity and 
density or velocity can be measured without knowing thick- 
ness. A pulse-echo terahertz velocity measurement is made 
by sending terahertz energy via a transceiver (device that has 
1 5 both a transmitter and a receiver) into and through a dielectric 
(insulating) material (such as the shuttle external tank thermal 
protection system sprayed on insulating foam) backed by a 
metallic (electrically -conducting) plate that reflects the tera- 
hertz energy back to the transceiver. The terahertz transceiver 
20 is separated from the dielectric sample by an air path. Velocity 
(V) values are calculated using the time delay (a) between the 
front surface echo (FS) and substrate reflection (BS) (with 
sample present). The novel pulse-echo method described 
herein for measuring velocity in a material sample uses echos 
25 off of the reflector plate without the sample present as well as 
the FS and BS echos with the sample present. 

With a dielectric sample present between the transceiver 
and the reflector plate, the pulse that travels from the trans- 
ceiver through the sample to the reflector plate (equivalent to 
30 the sample back surface position) and back to the transceiver 
is labeled BS and will be observed at time t\ Placing a dielec- 
tric sample in between the terahertz transceiver and the reflec- 
tor plate slows down the terahertz pulse as compared to its 
travel time in air. Thus, with the sample removed, the pulse 
35 that travels from the transceiver to the reflector plate and back 
to the transceiver is labeled M" and will be observed at an 
earlier time t M . For certain materials such as foams, by appro- 
priate manipulations and substitutions of equations, the 
acquisition of scans of the FS, BS, and M" (echo off reflector 
40 without sample present) echos, the conditioning of the FS 
echo thru amplification, DC subtraction, and software denot- 
ing, the fusing (combining) of FS, BS, and M M data sets 
(through use of software), and the subsequent calculation of 
time delays between echos, precision thickness -independent 
45 velocity images (that map micro structure) and micro struc- 
ture-independent thickness images (that map thickness) 
(through software) are obtained. 

The use of terahertz energy to simultaneously determine 
density and thickness variation in dielectric materials is new. 
50 The method is totally non-contact, very precise, and involves 
no fluid immersion. 

Prior to implementation of this method, it was not possible 
to separate out effects of thickness and microstructural varia- 
tion in time-of-flight images in totally non-contact, non-wa- 
55 ter-immersion fashion. No attempts to separate thickness and 
microstructural effects in terahertz time-of-flight images 
were noted in the literature. Ultrasonic methods to simulta- 
neously measure or characterize thickness and density (or 
variation as such) require water coupling. Additionally, the 
60 ultrasonic methods cannot be used for foam inspections due 
to the highly porous nature or highly cellular structure of 
foams. The terahertz method is totally non-contact, requires 
no coupling, and works in air. 

A process for non-destructive evaluation of a sample using 
65 a computer and a computer program implementing an algo- 
rithm is disclosed and claimed. The algorithm includes deter- 
mining the velocity of terahertz electromagnetic radiation in 



US 7,933,027 B1 


7 

the sample without prior knowledge of the thickness of the 
sample and determining the thickness of the sample without 
prior knowledge of the velocity of the terahertz electromag- 
netic radiation in the sample. The terahertz electromagnetic 
radiation is produced by a source spaced apart from the 
sample under test and propagated at the speed of light, c, in a 
medium located between the source and the sample. The steps 
include emitting and scanning terahertz electromagnetic 
radiation from the source. The terahertz radiation is applied to 
a plurality of scan points arranged according to x and y 
coordinates of the source. The plurality of scan points 
includes an area at least as large as the surface of the sample. 
Recording and storing in a substrate echo waveform data set 
file, on a scan point by scan point basis, the substrate echo 
(M M ) waveform of the terahertz electromagnetic radiation 
received from the substrate without the sample present is 
performed. The next step is placing the sample on the sub- 
strate. Next, the step of recording and storing in a back surface 
echo waveform data set file, on a scan point by scan point 
basis, the back surface echo (BS) waveform of the terahertz 
electromagnetic radiation received from the substrate with 
the sample present, f (t-prime, back surface echo (BS)). Then 
an additional step of recording and storing in a front surface 
echo (FS) waveform data set file, on a scan point by scan point 
basis, is undertaken. The front surface echo (FS) waveform of 
the terahertz electromagnetic radiation is received from the 
front surface of the sample. 

Once the aforementioned waveform files are generated, the 
step of recalling the stored substrate echo (M M ) waveform 
data set file, recalling the stored back surface echo (BS) 
waveform data set file, and recalling the stored front surface 
echo (FS) waveform data set file is performed. Next the step 
of opening and viewing, interactively, the substrate echo (M") 
waveform, the front surface echo (FS) waveform and the back 
surface (BS) echo waveform is performed. The substrate echo 
(M M ) waveform, the front surface echo (FS) waveform, and 
the back surface echo (BS) waveform for a single scan point, 
are superimposed in a display indicating amplitude and time 
base. The software next permits viewing, interactively, an 
algorithm-generated amplitude display of all scan points and 
associated substrate (M M ) waveforms, associated front sur- 
face (FS) waveforms and associated back surface echo (BS) 
waveforms, for selected scan points of the sample arranged 
according to x and y scan coordinates. The user may canvass 
any of the superimposed waveforms to obtain an idea or range 
of values for subsequent gating. Next, the software applies 
gating, interactively, the front surface echo (FS) waveform 
and the substrate echo (M M ) waveform based on ranges deter- 
mined from viewing selected scan points of the sample. The 
gating windows the FS and M M echoes of the scan points so 
that the windowed portions can be extracted and combined 
(fused) with the BS echoes at the same scan point. Condition- 
ing, if necessary of the front surface echo (FS) waveform is 
facilitated by the software. 

Fusing the stored substrate echo (M M ) waveform data set 
file, the stored back surface echo (BS) waveform data set file, 
and the stored front surface echo (FS) waveform data set file, 
together, producing a fused data file for each of the scan 
points, is performed by the software. Next the software facili- 
tates viewing, interactively, an algorithm-generated ampli- 
tude display of the fused waveforms associated with all scan 
points for selected scan points of the sample arranged accord- 
ing to x and y scan coordinates. 

Again, the software facilitates gating, interactively, the 
front surface (FS) echo waveform and gating, interactively, 
the back surface (BS) echo waveform of one of the fused 
waveform files based on ranges determined from viewing a 
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range of the selected fused waveforms. A calculation of the 
difference in time, 2r, between the front surface (FS) echo 
waveforms and the back surface (BS) echo waveforms on a 
scan point by scan point basis is made by the software and a 
5 2t file is generated and stored. 

Again, the software facilitates gating, interactively, the 
back surface (BS) echo waveform and gating, interactively, 
the substrate (M") echo waveform of one of the fused wave- 
form files based on ranges determined from viewing one of 
10 the selected fused waveforms. Next, a calculation of the dif- 
ference in time, At, between the front surface (BS) echo 
waveforms and the substrate (M M ) echo waveforms on a scan 
point by scan point basis is generated and stored as a At file. 
15 The generation of the 2 t and At files enables the step of 
determining and storing, on a scan point by scan point basis, 
the velocity of terahertz electromagnetic radiation in the 
sample without prior knowledge of the thickness of the 
sample. Further, generation of the 2 t and At files enables 
20 determining and storing, on a scan point by scan point basis, 
the thickness of the sample without prior knowledge of the 
velocity of the terahertz electromagnetic radiation in the 
sample. 

It is an object of the invention to simultaneously: (1) mea- 
25 sure the velocity of terahertz electromagnetic radiation in a 
material sample without prior knowledge of the thickness of 
the sample; and, (2) measure the thickness of a material 
sample using terahertz electromagnetic radiation in a material 
sample without prior knowledge of the velocity of the tera- 
30 hertz electromagnetic radiation in the sample. 

It is an object of the invention to provide a non-contact 
single-sided terahertz electromagnetic measurement and 
imaging method that simultaneously characterizes micro- 
structural (for example, spatially-lateral density) and thick- 
35 ness variation in dielectric (insulating) materials. 

It is an object of the present invention to provide a non- 
contact single-sided terahertz electromagnetic measurement 
and imaging method that simultaneously characterizes 
microstructural and thickness variation in dielectric (insulat- 
40 ing) materials. 

It is an object of the present invention to provide an inspec- 
tion method for current and future thermal protection systems 
and for other dielectric material inspection applications 
where microstructural and thickness variation require preci- 
45 sion mapping. 

It is an object of the present invention to provide an inspec- 
tion method which allows the separation of time-of-flight 
variations into its microstructural and thickness components. 

It is an object of the present invention to provide simulta- 
50 neous noncontact precision imaging of microstructural and 
thickness variation in dielectric materials using terahertz 
energy. 

It is an object of the present invention to provide simulta- 
neous noncontact precision imaging of microstructural and 
55 thickness variation in dielectric materials using terahertz 
energy using fused waveforms of terahertz energy from a 
fused data files produced by merging the FS, BS, and M" scan 
data sets for a set of foam blocks. 

It is an object of the present invention to provide an imag- 
60 ing of microstructural and thickness variation in dielectric 
materials using terahertz electromagnetic through the cre- 
ation of scan data arranged in an x-y matrix of scan points. 

It is an object of the present invention to provide a com- 
puter algorithm which enables mapping of scan data arranged 
65 in an x-y matrix which provides information about the veloc- 
ity/density of the sample under evaluation/test without know- 
ing the thickness of the sample under evaluation/test. 
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It is an object of the present invention to provide a com- 
puter algorithm which enables mapping of scan data arranged 
in an x-y matrix which provides information about the thick- 
ness of the sample under evaluation/test without knowing the 
density of the sample under evaluation. 5 

It is an object of the present invention to provide a system 
which includes a computer algorithm for creating back sur- 
face (BS) waveform data files, substrate echo (M M ) waveform 
data files, and front surface (FS) waveform data files. 

It is an object of the present invention to provide a system to 
which includes a computer algorithm which provides inter- 
active viewing of the back surface (BS) waveform data files, 
substrate echo (M M ) waveform data files, and front surface 
(FS) waveform data files for the determination of gating (win- 
dowing) the substrate echo (M M ) waveforms and the front 15 
surface (FS) waveforms and the subsequent fusing of the data 
files on a scan point by scan point basis for the creation of 
fused data files for each scan point. 

It is an object of the present invention to provide a system 
which includes a computer algorithm which provides gating 20 
the front surface (FS) echo and the back surface (BS) echo for 
the determination of a matrix of 2x data files on a scan point 
by scan point basis. 

It is an object of the present invention to provide a system 
which includes a computer algorithm which provides gating 25 
the substrate (M M ) echo and the back surface (BS) echo for the 
determination of a matrix of At data files on a scan point by 
scan point basis. 

It is an object of the present invention to provide a system 
which includes an image of the velocity of the terahertz 30 
radiation in the sample, and, hence the density of the sample 
without prior knowledge of the thickness of the sample. 

It is an object of the present invention to provide a system 
which includes an image of the thickness the sample without 
prior knowledge of the velocity of the radiation in the sample. 35 

It is an object of the present invention to provide a com- 
puter and computer software for recording and storing front 
surface (FS) echo files, substrate echo (M M ) files and back 
surface (BS) echo files and superimposing them in a time 
based file on a scan point by scan point basis. 40 

It is a further object of the present invention to gate (win- 
dow) superimposed data after inspecting the superimposed 
data and to produce a fused file for each scan point which 
includes a front surface (FS) echo waveform, a substrate (M M ) 
waveform, and a back surface (BS) echo waveform fused 45 
together. 

It is a further object of the present invention to gate a 
selected fused waveform for calculating 2 t and At on a scan 
point by scan point basis and then calculate and output veloc- 
ity (density) and thickness plots. 50 

These and other objects will be better understood when 
reference is made to the drawings, the description of the 
invention and claims which follow herein below. 

BRIEF DESCRIPTION OF THE DRAWINGS 55 

FIG. 1 is a schematic diagram of the reflection-mode tera- 
hertz methodology illustrating reflections received from the 
various interfaces and as an example of gating of the reflected 
signals. 60 

FIG. 2 is a plot of velocity versus density for foam. 

FIG. 2 A is a plot of velocity versus density for silicon 
nitride. 

FIG. 3 illustrates a schematic of the pulse-echo terahertz 
testing method and resulting waveforms, including BS (time 65 
to and from the back surface of the dielectric sample), FS 
(time to and from the front surface of dielectric sample), L 
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(distance between transceiver and sample), M M (pulse that 
travels from the transceiver to the reflector plate and back to 
the transceiver), d (sample thickness), f (t-prime) (travel time 
of the terahertz electromagnetic radiation to and from the 
substrate with the sample present, equal to 2 t l plus 2x), t" 
(t-double prime) (travel time of the terahertz electromagnetic 
radiation to and from the substrate without the sample 
present), At- (transmission time difference (f minus t") with 
the sample present and without the sample present), 2 t l 
(travel time of the terahertz electromagnetic radiation to and 
from the front surface of the sample) and 2x (time difference 
between the FS echo and the BS echo). 

FIG. 3 A illustrates a schematic similar to that illustrated in 
FIG. 3 with an additional thin dielectric material placed over 
the dielectric sample to create a dielectric mismatch. 

FIG. 4 is the uncertainty (in percentage) of the velocity 
(independent of the thickness) as a function of velocity, 2x 
and At. 

FIG. 5 is the uncertainty (in percentage) of the thickness 
(independent of the micro structural variations) as a function 
of thickness. 

FIG. 6 is a perspective view of a foam block sample set on 
an aluminum plate. 

FIG. 7 illustrates a composite image of scans of the foam 
block sample set of FIG. 6 indicating 2x values of the foam 
blocks according to the methodology of the invention. 

FIG. 7A illustrates an image of scans of the foam block 
sample set of FIG. 6 indicating thickness variations according 
to the methodology of the invention. 

FIG. 7B illustrates an image of scans of the foam block 
sample set of FIG. 6 indicating density variations according 
to the methodology of the invention. 

FIG. 8 illustrates a physically-measured density map in 
grams per cubic centimeter for a 6 by 15 set of foam blocks. 

FIG. 8A illustrates the density by shade of gray in grams 
per cubic centimeter for the physically-measured density map 
shown in FIG. 8. 

FIG. 9 illustrates a terahertz density map in grams per cubic 
centimeter for the 6 by 1 5 set of foam blocks of FIG. 8 derived 
from the velocity variations (determined independently of 
thickness) according to the methodology of the invention 
using the relationship between terahertz velocity and density 
for foam shown in FIG. 2. 

FIG. 9A illustrates the density by shade of gray in grams 
per cubic centimeter for the terahertz density map of FIG. 9. 

FIG. 10 illustrates a hand-measured thickness map in cen- 
timeters for the 6 by 15 set of foam blocks. 

FIG. 10A illustrates the thickness by shades of gray indi- 
cated in centimeters for the thickness map of FIG. 10. 

FIG. 11 illustrates a terahertz thickness map for the 6 by 1 5 
set of foam blocks (determined independently of velocity) 
according to the methodology of the invention. 

FIG. 11A illustrates the thickness by shade of gray in 
centimeters for the terahertz thickness map of FIG. 11. 

FIG. 12 is a schematic diagram of a process for measuring 
the velocity of terahertz electromagnetic radiation in a mate- 
rial sample without prior knowledge of the thickness of the 
sample, the terahertz electromagnetic radiation produced by a 
source (transceiver) spaced apart from the sample under test 
and propagated at the speed of light, c, in a medium located 
between the source (transceiver) and the sample, the sample 
residing on a substrate. 

FIG. 13 is a schematic diagram of a process for measuring 
the thickness of a material sample using terahertz electromag- 
netic radiation in a material sample without prior knowledge 
of the velocity of the terahertz electromagnetic radiation in 
the sample, the terahertz electromagnetic radiation produced 
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by a source (transceiver) spaced apart from the sample under 
test and propagated at the speed of light, c, in a medium 
located between the source (transceiver) and the sample, the 
sample residing on a substrate. 

FIG. 14 is a schematic of a process for simultaneously 
measuring the velocity of terahertz electromagnetic radiation 
in a material sample without prior knowledge of the thickness 
of the sample and for measuring the thickness of a material 
sample using terahertz electromagnetic radiation in a material 
sample without prior knowledge of the velocity of the tera- 
hertz electromagnetic radiation in the sample, the terahertz 
electromagnetic radiation produced by a source (transceiver) 
spaced apart from the sample under test and propagated at the 
speed of light, c, in a medium located between the source 
(transceiver) and the sample, the sample residing on a sub- 
strate. 

FIG. 15 is a typical “fused” waveform from a “fused” data 
file produced by merging the FS, BS, and M M scan data sets 
for a set of foam blocks. 

FIG. 16 is a schematic of a process (another example) for 
simultaneously measuring the velocity of terahertz electro- 
magnetic radiation in a material sample without prior knowl- 
edge of the thickness of the sample and for measuring the 
thickness of a material sample using terahertz electromag- 
netic radiation in a material sample without prior knowledge 
of the velocity of the terahertz electromagnetic radiation in 
the sample, the terahertz electromagnetic radiation produced 
by a source spaced apart from the sample under test and 
propagated at the speed of light, c, in a medium located 
between the source and the sample, the sample residing on a 
substrate. 

FIGS. 17-17C are a schematic of an exemplary computer 
implemented process for nondestructively determining tera- 
hertz electromagnetic radiation velocity without prior knowl- 
edge of thickness of a sample for determining thickness with- 
out prior knowledge of the terahertz electromagnetic 
radiation velocity in the sample. 

FIGS. 18-18A are another schematic of an exemplary com- 
puter implemented process for nondestructively determining 
terahertz electromagnetic radiation velocity without prior 
knowledge of thickness of a sample for determining thickness 
without prior knowledge of the terahertz electromagnetic 
radiation velocity in the sample. 

FIG. 19 is another schematic of an exemplary computer 
implemented process for nondestructively determining tera- 
hertz electromagnetic radiation velocity without prior knowl- 
edge of thickness of a sample for determining thickness with- 
out prior knowledge of the terahertz electromagnetic 
radiation velocity in the sample. 

FIG. 20 is an oscilloscope-like trace of a front surface (FS) 
echo waveform and a back surface (BS) echo waveform for 
the x=86, y=84 scan point coordinate. 

FIG. 20A is an oscilloscope-like trace of a substrate echo 
(M M ) waveform for the x=86, y=84 scan point coordinate. 

FIG. 20B is an oscilloscope-like trace of the front surface 
(FS) echo waveform and the substrate echo (M M ) waveform 
superimposed into the back surface (BS) echo waveform for 
the x=86, y=84 scan point coordinate. 

FIG. 20C is an image of the back surface (BS) echo wave- 
form illustrating the x-y scan points ranging from x=171 to 
y=171 allowing the user to canvas any of the scan points to 
produce superimposed files having the front surface (FS) 
echo waveform, the substrate echo (M") waveform superim- 
posed into and onto the back surface (BS) echo waveform. 

FIG. 20D is an image of the superimposed file for scan 
coordinates of x=20, y=l 50 and the corresponding oscillo- 
scope-like trace of the superimposed files having the front 
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surface (FS) echo waveform, the substrate echo (M M ) wave- 
form superimposed into and onto the back surface (BS) echo 
waveform. 

FIG. 20E is an image of the superimposed file for scan 
5 coordinates of x=86, y=84 and the corresponding oscillo- 
scope-like trace of the superimposed files having the front 
surface (FS) echo waveform, the substrate echo (M") wave- 
form superimposed into and onto the back surface (BS) echo 
waveform. 

to FIG. 20F is an image of the superimposed file for scan 
coordinates of x=150, y=20 and the corresponding oscillo- 
scope-like trace of the superimposed files having the front 
surface (FS) echo waveform, the substrate echo (M M ) wave- 
form superimposed into and onto the back surface (BS) echo 
15 waveform. 

FIG. 20G is a screen shot from the software program 
implementing signal condition of the front surface (FS) echo 
waveform and delaying, by a specified amount of time, the 
substrate echo (M H ) waveform to facilitate gating of the sub- 
20 strate echo (M M ) waveform. 

FIG. 20H is an oscilloscope-like trace of the time shifted 
substrate echo (M M ) waveform and gating of the front surface 
(FS) echo waveform and the time shifted substrate (M") 
waveform. 

25 FIG. 201 is a screen shot similar to FIG. 20G from the 
software program implementing signal condition of the front 
surface (FS) echo waveform and delaying, by a specified 
amount of time, the substrate echo (M") waveform to facili- 
tate gating of the substrate echo (M M ) waveform and in addi- 
30 tion indicating the “fuse gated waves” boolean button. 

FIG. 20J is an oscilloscope-like trace of the fused wave- 
form containing the front surface (FS) echo, the substrate 
echo (M M ) and the back surface (BS) echo. 

FIG. 2 OK is the oscilloscope-like trace of the fused wave- 
35 form containing the front surface (FS) echo, the substrate 
echo (M M ) and the back surface (BS) echo of FIG. 20J illus- 
trating gates placed over the front surface (FS) echo and the 
back surface (BS) echo. 

FIG. 20L is the oscilloscope-like trace of the fused wave- 
40 form containing the front surface (FS) echo, the substrate 
echo (M M ) and the back surface (BS) echo of FIG. 20J illus- 
trating gates placed over the substrate echo (M M ) and the back 
surface (BS) echo. 

FIG. 20M is a computer generated image of the thickness - 
45 independent velocity image generated by the computer pro- 
gram. 

FIG. 20N is a computer generated image of the density 
image corresponding to the velocity image of FIG. 20M gen- 
erated by the computer program. 

50 FIG. 21 is a screen shot of the program illustrating a his- 
togram of the density plot and a density image after the 
outliers have been removed. 

DESCRIPTION OF THE INVENTION 

55 

FIG. 1 is a schematic diagram 100 of the reflection-mode 
terahertz methodology illustrating reflections received off of 
the various interfaces and gating of the reflected signal. A 
transceiver 101 includes a femtosecond laser (70 MHZ rep. 
60 rate) 102 which generates short terahertz electromagnetic 
pulses and a receiver 103. Reflections will be received from 
the various interfaces 109, 104. Reflection from the metal 
substrate 104 will be the strongest. The horizontal dotted line 
from the echo shows a time gate 110 typically used during 
65 signal processing. The back surface 105 of the metal substrate 
is illustrated as is the beginning of the void 106 in the foam, 
silicon nitride or other dielectric 108. The void in this 
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example terminates 107 at the front surface of the metal 
substrate 104. The front surface of the foam, silicon nitride or 
other dielectric 109 is illustrated in FIG. 1 as is a graphical 
depiction of a gate 110 for signal analysis. 

Terahertz imaging is being used at NASA for nondestruc- 
tive evaluation of the Space Shuttle external tank thermal 
protection system sprayed-on foam insulation (SOFI). Gen- 
erally, the terahertz method is used in the pulse-echo c-scan 
configuration to map variations in the peak amplitude of the 
echo off of the metal substrate after it has traveled through the 
foam. Traditional c-scan imaging scales the peak amplitude 
values (to an 8- or 16-bit gray or color scale) at each scan 
location to form an image. 

An additional implementation of pulse-echo c-scan imag- 
ing involves mapping variations in the time-of-flight of a 
terahertz echo peak, or mapping the time delay between front 
surface and substrate (with the sample present) echos (FS, 
BS). This implementation concerns itself more with mapping 
thickness or global microstructural variation (such as physi- 
cal density variation) as opposed to discrete flaw detection. 
Time delay (with the sample present 2x) between the front 
surface echo (FS) and substrate echo (BS) is directly affected 
by thickness variation and terahertz velocity in the material 
according. See, FIG. 3. Here the designations (2 t) and (2d) 
(versus x and d) are used since the ultrasonic echo travels 
through the material thickness in the pulse-echo mode. 

Terahertz velocity is affected by variations in a volumetric 
microstructural property such as physical density as illus- 
trated in FIG. 2, similar to the way ultrasonic velocity 
responds to microstructural variation. 

Determining the relationship between velocity and density 
allows density maps to be obtained from velocity maps as set 
forth herein. Spatial variations in part thickness and/or spa- 
tially-lateral microstructural character will result in varia- 
tions in maps of 2 t. Analagous to a complex number having 
real and imaginary parts, 2 t images can be thought of as 
having thickness and microstructural components if both 
thickness and microstructural variation are present. 

A terahertz method which allows the separation of time- 
of-flight variations into its microstructural and thickness 
components is disclosed herein. This method is important 
because it determines the extent of microstructural variation 
in a part that also has thickness variation. Additionally, it 
provides a non-contact method for mapping thickness and/or 
density. FIG. 2 is a plot 200 of terahertz electromagnetic 
radiation velocity versus density for sprayed on foam. This 
enables the conversion of a given velocity into a respective 
density. FIG. 2A is a plot 200A of terahertz electromagnetic 
velocity versus density for silicon nitride. 

FIG. 3 illustrates a schematic 300 of the pulse-echo tera- 
hertz testing method and resulting waveforms (output volt- 
ages), including: BS (time to and from the back surface of 
dielectric sample); FS (time to and from the front surface of 
the dielectric sample); L (distance between transceiver and 
sample); M M (pulse that travels from the transceiver to the 
reflector plate and back to the transceiver), d (sample thick- 
ness), f (t-prime) (travel time of the terahertz electromagnetic 
radiation to and from the substrate with the sample present, 
equal to (2^) plus (2x); t" (t-double prime) (travel time of the 
terahertz electromagnetic radiation to and from the substrate 
without the sample present), At transmission time difference 
(f minus t") with the sample present and without the sample 
present; 2 1 1 (travel time of the terahertz electromagnetic 
radiation to and from the front surface of the sample); and, 2x 
(time difference between the FS echo and the BS echo). 

As stated previously, the terahertz method of inspecting 
metal reflector-backed dielectric materials is utilized to 
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simultaneously provide thickness-independent velocity (free 
of thickness effects) and microstructure-independent thick- 
ness (free of microstructure effects) images. A pulse-echo 
terahertz velocity measurement is made by sending terahertz 
5 energy via a transceiver (device that has both a transmitter and 
a receiver) into and through a dielectric (insulating) material 
backed by a plate (electrically -conducting, generally metal- 
lic) that reflects the terahertz energy back to the transceiver. 
The terahertz transceiver is separated from the dielectric 
to sample by an air path. 

The novel pulse-echo method described herein for measur- 
ing velocity in a material sample uses echoes off of the reflec- 
tor plate with (BS) and without the sample present (M"), as 
well as using the echo (FS) off of the sample front surface. 
15 The following steps illustrate how velocity (V) in a sample of 
thickness (d) is determined without prior knowledge of thick- 
ness. With a dielectric sample present between the transceiver 
and the reflector plate, the pulse that travels from the trans- 
ceiver through the sample to the reflector plate (equivalent to 
20 the sample back surface position) and back to the transceiver 
is labeled BS and will be observed at time t' where: t'=(2^ +2x) 

Referring to FIG. 3, the pulse-echo terahertz testing and 
resulting waveforms of FS and BS occur with the sample 
present. M" occurs without the sample present. (2 f A ) and (2x) 
25 are the pulse-echo time delays of the terahertz pulse from the 
transceiver to the sample front surface and from the sample 
front surface to the substrate with the sample present. 

Placing a dielectric sample in between the terahertz trans- 
ceiver and the reflector plate slows down the terahertz pulse as 
30 compared to its travel time in air. Thus, with the sample 
removed, the pulse that travels from the transceiver to the 
reflector plate (metal substrate) 104 and back to the trans- 
ceiver is labeled and will be observed at an earlier time t" 
where: 

35 

t ,, =(2t l +2d/c)) 

where, c, is the velocity of terahertz energy in air and, d, is the 
air gap equal to the sample thickness. 

The velocity of light at standard temperature and pressure 
40 was used for c in this investigation and is equal to 
0.02997055434 cm/psec. Subtracting the time t" measured 
without the sample from the time measured with the sample, 
t’, yields At, follows: 

45 At=t'-t"=(2x-2d/c) 

The thickness (d) of the sample can be determined in the 
pulse-echo configuration from: 

2d=(2x)V 

5° 

which is simply velocity times time through the sample in 
both the forward and reverse directions. Solving for “d” and 
rearranging yields an expression for the velocity: 

55 V=c(l-At/2x) 

As seen from the equation for velocity, sample thickness 
(d) is not a variable in the equation. Thus, this method does 
not require prior knowledge of sample thickness. If extended 
to multiple measurements across the sample (imaging), 
60 sample thickness variation effects are eliminated in the image 
allowing a true picture of microstructural variation for types 
of microstructural variation (such as density variation) that 
correlate with and will be revealed by velocity variation. For 
conventional time-of-flight imaging which does not separate 
65 velocity, V, and thickness, d, any thickness variation effects 
would corrupt the evaluation of microstructural variation (de- 
termined from velocity, V). Thus the new methodology 
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allows true characterization of microstructural variation (i.e., 
density variation) in a material structure that is also nonuni- 
formly thick. 

The derived equation, namely, V=c(l-At/2x), illustrates 
how the terahertz velocity in a dielectric material will be 5 
reduced fractionally from that in air by the factor: 

(l-At/2x). 

Further, rearranging At=f-t M =(2x-2<i/c), to solve for sample 
thickness, d, yields: 10 

d=c(2x-At)/2 

which allows the calculation of absolute material thickness 
without prior knowledge of velocity. If extended to multiple 
measurements across the sample (imaging), sample micro- 15 
structure variation effects are eliminated in the image allow- 
ing a true mapping of thickness variation. For conventional 
thickness mapping, microstructure variation effects would 
corrupt the evaluation of thickness variation. 

Thus, the new methodology allows true characterization of 20 
thickness variation in a material structure that is of nonuni- 
form microstructure. A key point of the methodologies dis- 
closed and claimed herein is that both thickness -independent 
velocity and micro structure-independent thickness images 
can be derived from the same set of scan information. 25 

In practice 2x is experimentally obtained from the pulse- 
echo time delay between the first front surface echo (FS) and 
substrate echo (BS) with the sample present. Either the time 
difference from FS peak location to BS peak location or 
cross -correlation of the waveforms of the two echoes can be 30 
used to obtain the 2x time delay. At is the pulse-echo time 
difference between the echos off the reflector plate with (BS) 
and without (M M ) the sample present, respectively. 

In fact, after the shuttle flight STS- 1 14, the ability to non- 
destructively detect crushed foam became a significant prior- 35 
ity. The microstructure-independent thickness mapping 
method can be used to identify and quantify areas of crushed 
(pushed-in) foam and precisely map thickness. The thick- 
ness-independent velocity method can be used to identify and 
quantify density variations in foam and other materials. It is 40 
worth noting that the previously -discussed ultrasonic meth- 
ods for thickness-independent velocity and micro structure- 
independent thickness require water coupling while no such 
coupling is needed for terahertz methods. The latter fact 
makes the terahertz method much more practical than the 45 
ultrasonic method for dielectric materials. 

The front surface echo (FS) from the dielectric material 
may be of very low signal-to-noise ratio (SNR) depending on 
the dielectric match between air and the sample. FIG. 3A 
illustrates a schematic similar to that illustrated in FIG. 3 with 50 
an additional thin dielectric material 301 placed over the 
dielectric sample to create a dielectric mismatch. If a good 
dielectric match exists, much of the terahertz energy will be 
transmitted into the sample and this presents somewhat of a 
problem. Additionally, the focal plane sensitivity of the tera- 55 
hertz method for samples of nonuniform thickness may result 
in the FS echo too far out of focus and thus reduce the signal 
to noise ratio (SNR) even further, thus limiting the thickness 
variation over which the method can be used. The approxi- 
mate time location of the front surface echo (FS) off the 60 
dielectric sample must be known “a priori” and the wavetrain 
examined manually by an oscilloscope to determine special 
post-processing needs for amplification and denoising the 
front surface echo (FS). For the space shuttle external foam, 
the FS echo can be as small as Viooth the amplitude of the BS 65 
(back surface of the sample) echo. This requires signal pro- 
cessing/conditioning steps of denoising and/or low -pass 
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(smoothing) filtering followed by amplification (software 
gain) at the time location(s) of the FS echo to clearly separate 
the FS echo from baseline noise. 

As stated above, to create a better dielectric mismatch 
situation in which more of the terahertz energy is reflected 
back to the receiving system while an ample amount is still 
transmitted into the sample, a sheet of very thin (250 jam) 
plastic transparency paper 301 can be placed onto the sample. 
See, FIG. 3A. This method can be used to locate the front 
surface echo (FS) locations prior to scanning, or in situations 
where it can be tolerated during actual scanning, will provide 
front surface echos having much greater signal to noise ratios. 
Also, knowledge of the distance between scanner head and 
sample top surface, velocity of terahertz in air (speed of light), 
and any post- or pre-trigger delays should allow calculation of 
approximate front surface echo time location(s). 

FIG. 4 is a graph 400 of the uncertainty (in percentage) of 
the velocity (independent of the thickness) as a function of V 
(terahertz velocity in the material), 2x and At. The precision 
(uncertainty) in the thickness-independent velocity due to the 
random errors in the measurements of the variables At and 2x 
was determined by the above equations and standard variance 
relation. Uncertainty in c was ignored and using typical val- 
ues of At«6 gsec, 2x=200 psec, Sampling Rate=6.4 THz, 
c=0.02997055434 cm/psec, and V«0.0290 cm/psec, gives 
Uv^O.Ol percent (uncertainty of the velocity in percent). FIG. 
4 illustrates uncertainties (in percent) for three velocities, V, 
0.01 cm/psec; 0.03 cm/psec; and, 0.05 cm/psec. 

FIG. 5 is a graph 500 of the uncertainty (in percentage) of 
the thickness (independent of the microstructural variations) 
as a function of thickness, d. Similarly, the precision (uncer- 
tainty), Ud, of the thickness measurement is a function of the 
thickness as illustrated in FIG. 5. For the foam samples stud- 
ied, and using typical values for SR=6.4 THz, 
c=0.02997055334 cm/psec, and thicknesses of approxi- 
mately 3 to 5 cm, Ud=0.035 to 0.025 percent (uncertainty of 
the thickness in percent). 

FIG. 6 is a perspective view 600 of a foam block step wedge 
sample set on an aluminum plate 601. The first foam block 
602, the second foam block 603, the third foam block 604, 
and the fourth foam block 605 are illustrated and arranged 
with increasing thickness and density from left to right. The 
foam block configuration of FIG. 6 is ordered as step wedges 
such that density and thickness variation results in an additive 
effect on 2x. 

FIG. 7 illustrates a composite image 700 of scans of the 
foam block sample set of FIG. 6 indicating 2x values of the 
foam blocks according to the methodology of the invention as 
a function of thickness and density. FIG. 7 shows 2x, FIG. 7 A 
shows microstructure-independent thickness image and FIG. 
7B shows thickness-independent density image for the foam 
blocks of FIG. 6, using the novel terahertz method described 
in this invention. The density image was derived from thick- 
ness-independent velocity images using the relationship 
between terahertz velocity and density for foam shown in 
FIG. 2. This sample has well-defined thickness and density 
variation that have additive effects on the 2x image. It can be 
seen that the method is able to separate the thickness and 
density variation components as illustrated in FIGS. 7 A and 
7B. 

FIG. 7 A illustrates an image 700 A of scans of the foam 
block sample set of FIG. 6 indicating thickness variations 
according to the methodology of the invention. FIG. 7B illus- 
trates an image 70 0B of scans of the foam block sample set of 
FIG. 6 indicating density variations according to the method- 
ology of the invention. The method of this invention separates 
thickness and density effects on time delay (2x) between FS 
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and BS echos. Increasing thickness and increasing density 
from right-to-left in the blocks of FIG. 6 provides an additive 
effect in terms of increasing 2x from right-to-left. Scan and/or 
analog-to-digital conversion jitter (zigzag gray level pattern) 
is apparent in the density image of FIG. 7B. 

The terahertz experimental setup used in connection with 
the foam blocks of FIG. 6 was a broadband 1 THz scan system 
and included the following experimental and signal process- 
ing parameters: 

Focus: (At substrate, or 3 cm above substrate for wedge 
samples) 

Typical received bandwidth points (THz) (Full Width Half 
MaxfyO.l to 0.3 

Data Acquisition Rate (THz): 6.4 

Waveform Length Acquired (p sec/points): 320 psec/2048 
points 

Waveform Acquisition Rate (scanpoints/sec): «10 

Collinear source-detector: Yes 

Spatial Resolution (at Full Width Half Max of Point Spread 
Function): 0.5 cm 

Signal Acquisition Width of dynamic (peak-centered) 
gates for time delay computations (psec): 25 to 100 psec 

Scan Increment (cm): 0.2 

Samples were placed on an aluminum plate as shown in 
FIG. 6 and scanned with scan increments in the X and Y 
direction of 0.2 cm. The minimum number of scans required 
to obtain thickness -independent velocity and micro structure- 
independent thickness is two. One scan obtains FS and BS 
echos (with sample present) and the second scan obtains the 
M" echo (without sample present). Separate scans for FS and 
BS can be performed if sample thickness is too large to allow 
simultaneous capture of both of the echos in the 320 psec/ 
2048 point window using the 6.4 THz sampling rate. Only 
two scans were required when the thicknesses were <4 cm 
and FS and BS could be captured in one scan. For each sample 
set, the scans were then fused (combined using software) such 
that FS, BS and M M echos were placed in a single wavetrain of 
640 psec/4096 points length with time relationships between 
the echos preserved. This occurs at each scan location to 
create the new fused data set. Precise time delays 2x and At 
were determined using cross-correlation between the echos’ 
waveforms. Phase relationships were examined for: (1) FS 
compared to BS; and, (2) BS compared to M M . All waves 
appeared to be in phase for the analysis made. If echos are 
in-phase with respect to each other, the time occurrence of the 
maximum in the correlation function was used to calculate 
time delay. If the echos were phase-inverted, the time occur- 
rence of the minimum in the correlation function should be 
used to calculate time delays. 

A typical “fused” waveform from a “fused” data file is 
illustrated in FIG. 15. FIG. 15 is a typical “fused” waveform 
1500 from a “fused” data file produced by merging the FS, 
BS, and M M scan data sets for a set of foam blocks. 

The fused data file is produced by merging the FS, BS, and 
M M scan data sets for foam block sets. FS and BS occur with 
the sample present. M M occurs without the sample present. 
For visualization purposes, M" has been artificially shifted to 
the left an additional 40 psec in FIG. 15 to avoid overlap 
between the echoes. FS has been denoised, amplified, and a 
DC component has been subtracted therefrom so as to allow 
2x time delay calculation. A 25 to 1 00 psec gate (window) was 
applied to account for variations in FS echo position due to 
thickness variations in the sample. The gated region contain- 
ing the FS echo was denoised using a wavelet process, then 
amplified by 10 to 40x, followed by subtraction of the DC 
component. The denoising process used the debauchies 05 
mother wavelet principle. The resultant FS echo was quite 
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useable as shown in FIG. 15. Two Tau (2x) and At are deter- 
mined using the fusing process and the entire waveforms are 
either cross -correlated or the peaks are precisely identified 
enabling the measurement of the time between them. A com- 
5 puter process with an appropriate algorithm is used to calcu- 
late 2x and At. 

FIG. 16 is a schematic 1600 of another process for simul- 
taneously measuring the velocity of terahertz electromag- 
netic radiation in a material sample without prior knowledge 
1 0 of the thickness of the sample and for measuring the thickness 
of a material sample using terahertz electromagnetic radia- 
tion in a material sample without prior knowledge of the 
velocity of the terahertz electromagnetic radiation in the 
1 5 sample, the terahertz electromagnetic radiation produced by a 
source spaced apart from the sample under test and propa- 
gated at the speed of light, c, in a medium located between the 
source and the sample. The process includes the steps of: 
merging the FS, BS, and M M scan data sets into a fused data 
20 file 1601; producing a “fused” waveform 1602; calculating 
time delays between echos 1603; determining precision 
thickness-independent velocity images that map micro struc- 
ture 1604; and, determining precision micro structure-inde- 
pendent thickness images that map thickness 1605. The steps 
25 of calculating the time delays and determining precision 
thickness and micro structure are performed by the step of 
cross-correlating the entire waveform of FS, BS and M" 1606. 
Alternatively, the steps of calculating the time delays and 
determining precision thickness and micro structure are per- 
30 formed by precisely identifying peaks of FS, BS and M" 
signals 1607. 

If a poor dielectric mi smatch occurs between the air and the 
sample, the step of gating and processing the FS signal prior 
to merging the data into a fused data file 1608 is performed. If 
35 the signal to noise ratio is low then the step of amplifying the 
processed FS signal prior to merging the data into a fused data 
file 1609 is performed. Further the step of applying a dielec- 
tric sheet on the FS, followed by gating and processing the FS 
signal prior to merging the data into a fused data file 1610 may 
40 optionally be performed to enhance the signal processing. 

FIG. 8 illustrates 800 a physically -measured density map 
in grams per cubic centimeters for a 6 by 1 5 set of foam blocks 
with non regular thickness and density. The set of 6 by 15 
foam blocks had dimensions of about 5 by 5 by 5 cm, with 
45 minor but non regular thickness variation (±0.1 cm). The 
blocks were of various densities ranging from about 0.042 to 
0.054 g/cm 3 “(on the order of 20 percent)” measured from 
mass and dimensional measurements and were arranged ran- 
domly. FIG. 8 A illustrates 8 00 A the density by shade of gray 
50 in grams per cubic centimeter for the physically-measured 
density map shown in FIG. 8. FIG. 9 illustrates a terahertz 
density map 900 for the same 6 by 15 set of foam blocks 
derived from the velocity variations (determined indepen- 
dently of thickness) according to the methodology of the 
55 invention using the relationship between terahertz velocity 
and density for foam shown in FIG. 2. FIG. 9 A illustrates 
9 00 A the density by shade of gray in grams per cubic centi- 
meter for the derived terahertz density map shown in FIG. 9. 

Reference numeral 801 represents an area of the physically 
60 measured density map and reference numeral 901 represents 
an area of the mapped terahertz density plot for the same set 
of foam blocks. Reference numeral 802 represents an area of 
the physically-measured density shown in FIG. 8 of the set of 
foam blocks and reference numeral 902 represents an area of 
65 the mapped terahertz density plot for the same set of foam 
blocks. Viewing lighter and darker areas (801, 901, 802, 902) 
in the images of FIGS. 8 and 9, it is clear that the physically- 
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measured density variation agrees quite closely with that 
derived from the thickness -independent velocity. 

FIG. 10 illustrates a hand-measured thickness map 1000 in 
centimeters for the same 6 by 15 set of foam blocks. FIG. 11 
illustrates a terahertz thickness image 1100 for the 6 by 15 set 
of foam blocks (determined independently of velocity) 
according to the methodology of the invention. FIG. 10A 
illustrates the thickness in centimeters by shade of gray in a 
bar graph 1000A. FIG. 11 A illustrates the thickness 1100A 
by shade of gray in centimeters for the terahertz thickness 
map of FIG. 11. As stated previously the blocks vary in 
thickness from 5 cm by ±0.1 cm. Reference numeral 1001 is 
an ellipse indicating an area of the hand measured thickness 
map to be compared to a terahertz thickness map and refer- 
ence numeral 1101 is an ellipse indicating an area of the 
terahertz thickness map to be compared to a physically mea- 
sured thickness map. Note the excellent correlation between 
dark and light areas in both images. Ellipses 1001 and 1101 
denote the identical area under examination. Good correla- 
tion within ellipses 1001 and 1101 between light and dark 
areas is observable. Dark scatter spots in the terahertz gener- 
ated thickness are due to the presence of an additional echo 
within the signal processing gate that results in improper 
cross-correlation delay calculation. This additional echo is 
likely due to the presence of extra material on the surface. 
These scatter spots are also in the thickness -independent 
velocity image of FIG. 9 but blend in better as they cause 
variation in that image in the same “direction” as actual 
velocity variations. 

FIG. 12 is a schematic diagram 1200 of a process for 
measuring the velocity of terahertz electromagnetic radiation 
in a material sample without prior knowledge of the thickness 
of the sample, the terahertz electromagnetic radiation pro- 
duced by a source (transceiver) spaced apart from the sample 
under test and propagated at the speed of light, c, in a medium 
located between the source and the sample. The steps of the 
process may be performed in any desired or necessary order. 
The process includes the steps of: emitting pulsed terahertz 
electromagnetic radiation from the source 1201 (the radiation 
may be pulsed or it may be continuous); measuring the travel 
time of the terahertz electromagnetic radiation to and from 
the substrate without the sample present, t" (t-double prime), 
1202; placing the sample on the substrate 1203; measuring 
the travel time of the terahertz electromagnetic radiation to 
and from the substrate with the sample/dielectric present, t' 
(t-prime) 1204; subtracting t" from t' to determine the trans- 
mission time difference with the sample present and without 
the sample present, At, 1205; measuring the travel time of the 
terahertz electromagnetic radiation to and from the sample, 
2t l , 1206; determining the time difference between the FS 
echo and the BS echo, 2x, 1207; placing a sheet of plastic 
transparent paper onto the sample to create an adequate 
dielectric mismatch 1208; dividing At by 2 t and determining 
the quotient 1209; subtracting the quotient from 1 to obtain a 
factor 1210; multiplying the factor by, c, to determine the 
velocity of the terahertz electromagnetic radiation in the 
sample 1211 ; determining the micro structural variation of the 
sample 1212; determining the density of the material 1213; 
evaluating, in a plurality of locations, the sample for micro- 
structural variations; and mapping the variations by location 
1214; controlling the spacing between the source and the 
substrate 1215; and, determining, initially, the approximate 
time location of the sample from the terahertz radiation 
source 1216. Additionally, the step of determining the micro- 
structual variation of the sample may include the determina- 
tion of a parameter other than density. In the case of a poor 
dielectric mismatch between the air and the sample to be 
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inspected, the process can include a step of placing a sheet of 
plastic transparent paper onto the sample to create an 
adequate dielectric mismatch. 

FIG. 13 is a schematic diagram 1300 of a process for 
5 measuring the thickness of a material sample using terahertz 
electromagnetic radiation in a material sample without prior 
knowledge of the velocity of the terahertz electromagnetic 
radiation in the sample, the terahertz electromagnetic radia- 
tion produced by a source spaced apart from the sample under 
to test and propagated at the speed of light, c, in a medium 
located between the source and the sample. The process for 
measuring the thickness using terahertz electromagnetic 
radiation includes: emitting terahertz electromagnetic radia- 
tion from the source 1301 (the radiation may be pulsed or it 
1 5 may be continuous) ; measuring the travel time of the terahertz 
electromagnetic radiation to and from the substrate without 
the sample present, t" (t-prime), 1302; placing the sample on 
the substrate 1303; measuring the travel time of the terahertz 
electromagnetic radiation to and from the substrate with the 
20 sample present, t' (t-prime), 1304; subtracting t" from t' to 
determine the transmission time difference with the sample 
present and without the sample present, At, 1305; measuring 
the travel time of the terahertz electromagnetic radiation to 
and from the sample, 2 t l9 1306; determining the time differ- 
25 ence between the FS echo and the BS echo, 2x, 1307; sub- 
tracting At from 2x to obtain a subtraction result 1308; divid- 
ing the subtraction result by 2 to obtain a quotient 1309; 
multiplying the quotient by, c, to obtain the thickness of the 
sample 1310; and, evaluating, in a plurality of locations, the 
30 sample for thickness variations; and mapping the thickness 
variations by location 1312. 

FIG. 14 is a schematic 1400 of a process for simultaneously 
measuring the velocity of terahertz electromagnetic radiation 
in a material sample without prior knowledge of the thickness 
35 of the sample and for measuring the thickness of a material 
sample using terahertz electromagnetic radiation in a material 
sample without prior knowledge of the velocity of the tera- 
hertz electromagnetic radiation in the sample, the terahertz 
electromagnetic radiation produced by a source spaced apart 
40 from the sample under test and propagated at the speed of 
light, c, in a medium located between the source and the 
sample. The process for simultaneously measuring the veloc- 
ity of terahertz electromagnetic radiation in a material sample 
without prior knowledge of the thickness of the sample and 
45 for measuring the thickness of a material sample using tera- 
hertz electromagnetic radiation in a material sample without 
prior knowledge of the velocity of the terahertz electromag- 
netic radiation in the sample includes the steps of: emitting 
terahertz electromagnetic radiation from the source 1401 (the 
50 radiation may be pulsed or it may be continuous) ; measuring 
the travel time of the terahertz electromagnetic radiation to 
and from the substrate without the sample present, t" 
(t-double prime) 1402; placing the sample on the substrate 
1403; measuring the travel time of the terahertz electromag- 
55 netic radiation to and from the substrate with the sample 
present, t' (t-prime) 1404: subtracting t" from f to determine 
the transmission time difference with the sample present and 
without the sample present, At, 1405; measuring the travel 
time of the terahertz electromagnetic radiation to and from 
60 the sample, 2t l lain, 1406; determining the time difference 
between the FS echo and the BS echo, 2x, 1407; dividing At 
by 2x and determining the quotient 1408; subtracting the 
quotient from 1 to obtain a factor 1409; multiplying the factor 
by, c, to determine the velocity of the terahertz electromag- 
65 netic radiation in the sample 1410; subtracting At from 2x to 
obtain a subtraction result 1412; dividing the subtraction 
result by 2 to obtain a quotient 1413 ; multiplying the quotient 
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by, c, to obtain the thickness of the sample 1414; and, evalu- 
ating, in a plurality of locations, the sample for microstruc- 
tural variations and for thickness variations, and mapping the 
microstructural and thickness variations by location 1416. 

Operation of the computer interface (NDE Wave And 
Image Processor) to create fused files, determine and store 
data files, and output velocity/density images and thickness 
images is as follows: 

Creating a Fused File 

A computer program for use with a computer and computer 
display has be created using Lab VIEW® software. First, the 
user must select a stored waveform file to open from an 
opening dialog box in the NDE Wave and Image Process 
program opening screen. Open a back surface file when 
prompted to open waveform file. The underlying code of this 
interface uses Lab VIEW® functions for the interface, code 
for dialog box, and basic Lab VIEW® functions for opening a 
file, and parsing information from binary files. Lab VIEW® is 
a registered trademark of National Instruments. Upon open- 
ing a file, a peak-to-peak amplitude image is created and 
shown in an image container and the waveform associated 
with the center point of the scan is shown in the time-domain 
graph. FIG. 20 illustrates a front surface (FS) echo waveform 
and a back surface echo waveform. FIG. 20D illustrates 
image windows for superimposed files. 

Next, in the menu, click on density/thickness study, under 
the miscellaneous menu. After “density/thickness study” is 
selected, controls will appear on the main window in the area 
where the image container was and the image container box 
disappears. The user is then prompted to open the three files 
containing echoes to be fused, the first of which is back 
surface file which was previously opened. The user then des- 
ignates the substrate echo (M M ) and front surface (FS) echo 
files next, and then the user selects “ok” to superimpose the 
files on the same oscilloscope-like trace. The back surface 
(BS) echo file and front surface (FS) echo files might actually 
be the same files but both files have to be entered into the 
respective dialog boxes when prompted to do so. The pro- 
gram prompts the user to enter a new file name. 

This file will consist of the three echoes fused together and 
needed for thickness -independent velocity and microstruc- 
ture-independent thickness. A new file is created and three 
files are superimposed in an oscilloscope trace and might look 
something like that shown by way of example in FIG. 2 OB. 
The echo files are not yet fused into a new file at this point but 
are superimposed on the same trace as the waveforms in the 
three files are recalled and positioned appropriately for visu- 
alization. See, FIG. 20B. 

The following steps result in actual fusing of the echoes 
into one file. The software prompts the user to interactively 
select “gate” and then select “multigate” under gate type 
menu. Adjust substrate echo offset delay as needed to space it 
away from the back surface echo. FIG. 20G illustrates an 
interactive selection of the time shift for the substrate echo 
(M M ) waveforms to place it away from the back surface (BS) 
echo wave prior to gating and fusing the waveforms together. 
Substrate echo will have higher amplitude than back surface 
echo. Generally, add (-) delay so as to move the substrate 
echo to the left away from the back surface echo. Referring to 
FIG. 20H, the user places gate 1, 2013, centered over the 
substrate echo (M M ) waveform 2013, and then mouse clicks 
the “Add Gate” button. The user then places gate 2, 2012, to 
the time location where front surface (FS) echo waveform is; 
apply digital filtering, wavelet denoising, and wavelet recon- 
struction as needed, adjust FS gate amplification factor, 
Z-axis move distance if refocusing on FS occurred (+value), 
as needed. See FIGS. 2000H and I. Once gate 2, 2012, is 
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placed the user mouse clicks “Add Gate” button. See FIGS. 
20G and 201 wherein the user of the software is prompted to 
select boolean buttons for signal conditioning the front sur- 
face (FS) echo waveform/signal. The front surface (FS) echo 
5 waveform may be in need of conditioning as it is sometimes 
a weak signal if the dielectric mismatch is not significant. 

Modify c (cm/sec) depending on whether experimental 
conditions are for electromagnetic waves in air or ultrasound 
in water. If the front surface (FS) echo waveform is very weak, 
10 first you must know approximate time location of FS echo in 
FS file since amplitude is almost in the noise and you need to 
properly gate. 

The user then mouse clicks “Fuse Gated Waves,” 2015, 
1 5 FIG. 201, which will place substrate echo (M M ) waveform and 
front surface echo (FS) waveform into the back surface (BS) 
echo waveform or wavetrain that has been expanded to 
double the size of the original wavetrain (egs. original=2048 
(320 psec), fused=4096 points (640 psec)). The code for 
20 fusing involves opening the three files and assembling the 
echoes using proper timing information and expanding the 
total number of positions in the waveform to twice the length 
of the individual waveform. 

Instruction for using the fused files follows. 

25 The software produces an image processing window when 
the file fusing procedure is completed and it will contain a 
peak-to-peak image. The cursor may be moved across the 
image processing window to recall the superimposed waves 
as needed. In effect, the user is canvassing areas of the sample 
30 to view the fused files for information as the range for later 
gating of the files. The time arrangement between substrate 
echo wave form (M M ) and back surface (BS) echo waveform 
will depend on whether evaluation involves electro -magnetic 
waves in air or ultrasound in water. For waves in air, the 
35 substrate echo M" waveform is always expected to occur 
earlier in time as the speed of light in air will be faster than that 
through a material. 

The next step includes opening a fused echo waveforms file 
containing fused waveforms for each scan point. Using the 
40 mouse, the user clicks on the image container window to 
spawn the image processing window. The waveforms in the 
time-domain graph in the main window should now have 
three distinct echoes based on the original gating of the front 
surface echo (FS) waveform and the substrate echo (M M ) 
45 waveform to form the fused file. An oscilloscope trace is 
displayed and should look something like FIG. 20J after 
fusion. The user can move and release a cursor in the image 
processing window to recall fused waveforms at any scan 
point. 

50 Gates are now used again and the selection “use gate” 
should be made as should the selection “multi gate” under gate 
type menu. Two gates should be placed as follows: gatel, 
2017, is placed on the front surface (FS) waveform echo and 
gate 2, 2018, is placed on the back surface (BS) waveform 
55 echo. See FIG. 20K. Form the delay image, an image of 2 t, 
between gates 1 and 2 (“peak to peak”, or “cross correlation”) 
under Form Image of Type/Calculate Time Domain Velocity- 
Thickness menus. If multiple front surface (FS) waveform 
times were gated during file fuse procedure, the front surface 
60 (FS) waveform gate (gate 1 , 2017) should be made as long as 
needed to cover the movement of the front surface (FS) wave- 
form echo determined by canvassing the waveforms in FIGS. 
20D-20F. Then the image should be saved as a .bin or .txt file. 
Saving as .txt or .bin file is done using basic Lab VIEW 
65 functions, “.bin” files are binary files and “.txt” files are text 
files and can be opened by spreadsheet programs to view 
outside of the Lab View program. 
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Referring to FIG. 20L, next form the At delay image 
between back surface echo (BS) waveform and substrate echo 
(M M , M M ) waveform M M echoes. First select using a mouse the 
choice “multi gate” under “Gate Type” menu to reset the gates 
selected. The gates should now be placed as follows: gatel 5 
(2019) on substrate echo (M M , M M ) waveform echo and gate 2 
(2020) on back surface (BS) echo waveform. Form the delay 
image (using either the peak to peak or cross -correlation 
methods), then save the image as a .bin or .txt file. 

To form thickness-independent velocity or microstructure 1C 
independent thickness image perform the following steps. 

Open the fused file and Select either Thickness -Indepen- 
dent Velocity or Micro structure-independent Thickness 
Image Type under Form Image of Type Calculate Velocity 
Thickness. Open either .bin or .txt files for 2x and At image 15 
files. The software algorithm uses the following equations 
and solves for V, velocity and d, thickness, as follows: 

Thickness-independent velocity is given by the following 
equation: 

c(2r- At) 
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Micro structure-independent Thickness is given by the fol- 
lowing equation: 


FIG. 20M is a velocity image 2000M which represents the 
thickness-independent velocity in gray shade according to the 
x,y coordinates selected and is then compared to the bar graph 35 
(gray scale ramp) 2021. Reference numeral 2022 illustrates 
the present x-y coordinate in the approximate center of FIG. 

20. To convert the images to a density image, the software 
prompts selection/conversion 2025 as illustrated in FIG. 
20M. FIG. 20N illustrates 2000N the density corresponding 40 
to the x,y coordinates and reference numeral 2024 illustrates 
the present x, y coordinate. 

To obtain the density maps/images from thickness -inde- 
pendent velocity maps/images, select the box on the graphical 
user interface “Convert to Density” which spawns a dialog 45 
prompting the entry of a slope (m) and intercept (b) of the best 
linear fit for the velocity -physical density relationship that 
has been determined through experiments and analysis of the 
material under test/evaluation. 

The final density image may require post-processing using 50 
another software module within NDE Wave & Image Proces- 
sor if background was obtained in the scan or if outlier data 
points from bad data are obtained. For example, contrast 
expansion/outlier removal of the sample area may need to be 
performed using histogram range-based contrast expansion 55 
module. Quantitative density and thickness maps are 
obtained. See FIG. 21. 

One exemplary process for non-destructive evaluation of a 
sample using a computer, a computer display and a computer 
program includes implementing an algorithm for determin- 60 
ing the velocity of terahertz electromagnetic radiation in the 
sample without prior knowledge of the thickness of the 
sample and determining the thickness of the sample without 
prior knowledge of the velocity of the terahertz electromag- 
netic radiation in the sample. The terahertz electromagnetic 65 
radiation is produced by a source spaced apart from said 
sample under test and propagated at the speed of light, c, in a 


24 

medium located between the source and the sample, and is set 
forth in an example schematically illustrated in FIG. 18. 

The first exemplary process comprises the steps of: emit- 
ting and scanning terahertz electromagnetic radiation from 
the source, the terahertz radiation applied to a plurality of 
scan points arranged according to x and y coordinates of the 
scanning, the plurality of scan points includes an area at least 
as large as the surface of the sample, 1801. The scan points 
may be arranged in an x-y matrix of 171 data points by 171 
data points (or any arrangement of data points as dictated by 
the sample size and scan resolution desired). 

This example of the process includes recording and stor- 
ing, in a substrate echo wavefonn data set file, on a scan point 
by scan point basis, a substrate echo (M M ) waveform of the 
terahertz electromagnetic radiation received from the sub- 
strate without the sample present, and recording and storing 
the front surface (FS) waveform and the back surface (BS) 
waveform of the terahertz electromagnetic radiation received 
with the sample present, 1802. So, for example, if there are 
171 2 data points then there are 29,241 data points and each 
data point or data area has associated therewith a substrate 
echo (M M ) waveform, a back surface (BS) echo wave form 
and a front surface (FS) echo waveform. These waveforms are 
superimposed on a time line (abscissa) and display amplitude 
with respect to time. See FIG. 20B for an illustration of 
superimposed waveforms for scan point x=86, y=84. Refer- 
ence numeral 2001 is the back surface (BS) echo waveform 
and the front surface echo waveform is indicated by reference 
numeral 2002 in FIGS. 20, 20A and 2 0B The substrate echo 
(M M ) waveform is illustrated in FIGS. 20A and 20B. There- 
fore, the step of superimposing the substrate echo (M M ) wave- 
form, the front surface (FS) echo waveform, and the back 
surface (BS) echo waveform, 1803 accomplishes the super- 
position. Gating (2013, 2012) the substrate echo (M M ), and 
gating and conditioning the front surface (FS) echo waveform 
are performed in the program as a precursor to fusing the echo 
waveforms together. See FIG. 2 OH which illustrates placing 
the first gate 2013 and the second gate 2012 over the substrate 
echo (M M ) waveform and the front surface (FS) echo wave- 
form, respectively. 

After superposition (FIG. 2 0B) has been performed an 
image container (FIG. 20C, display on the graphical user 
interface) appears and the user may interactively with the 
computer mouse or a touch screen interface canvas the sample 
under evaluation to determine the characteristics of the wave- 
forms on a point by point basis (x,y coordinates). This can- 
vassing or inspection gives the software user an idea of the 
breadth or variance of the front surface (FS) echo waveform 
over a given range. This range then determines the width of 
the gate over the front surface (FS) echo waveform. The range 
of the superimposed waveforms are illustrated in FIGS. 20D, 
20E, and 20F. Movement of the superimposed waves is illus- 
trated in comparing FIGS. 20D, 20E and 20F. FIG. 20D 
corresponds to x=20, y=150, FIG. 20E corresponds to x=86, 
y=84, and FIG. 20F corresponds to x=150, y=20. A single 
waveform is illustrated as an oscilloscope trace and it is this 
single waveform that is gated. See FIG. 20G which illustrates 
delaying the substrate echo (M M ) waveforms 2003 as indi- 
cated by the time shifted substrate echo (M M ) waveform 
2003A. The substrate echo (M M ) waveform is also gated. 

Signal conditioning of the front surface echo (FS), namely, 
applying digital filtering, wavelet denoising, and/or wavelet 
reconstruction, adjusting the front surface (FS) gate amplifi- 
cation factor, Z-axis move distance if refocusing on FS 
occurred (+value), FS Subtraction Offset, as needed may be 
applied. See FIG. 201, reference numeral 2014 indicating 
possible selections of signal conditioning processes. 
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The step of fusing the substrate echo (M M ) waveform and 
the front surface echo (FS) waveform into the back surface 
(BS) echo wavefonn is performed 1804 and selecting signal 
conditioning is performed on the front surface echo (FS) 
waveforms for each of the data points as it is fused with its 
respective substrate echo (M M ) waveform and the back sur- 
face (BS) waveforms. Fusing is accomplished after gating the 
waveforms by selecting the “fuse gated waves” action, 2015. 
See FIG. 201. Once all of the echo waveforms are fused into 
a fused data set the next step of this example is to gate the echo 
waveforms and determine the time difference between the 
front surface echo (FS) waveform and the back surface (BS) 
echo waveform, 2x, and storing values of 2 t for each scan 
point, 1805. The gates are interactively applied with the user 
of the software program, computer and computer display, 
actively dragging, dropping and adjusting the width (in pico 
seconds) of the gate 2017 over the respective front surface 
(FS) echo waveform and the gate 2018 over the back surface 
(BS) echo waveforms. Prior to selecting the position and 
width of the gates of the front surface (FS) waveform and the 
back surface waveform (BS) waveforms, the user may canvas 
the fused file (FIG. 20J) over the extent of the sample under 
evaluation to assist the user in deciding how wide the gates 
should be in time. The canvassing of the fused files is done in 
the image container in windows similar to those illustrated in 
FIGS. 20D-F. 

Next, the step of this example includes interactively gating 
the echo (2019, 2020) waveforms and determining the time 
difference between the substrate echo wave (M M ) waveform 
and the back surface (BS) waveform, At, and storing values of 
At for each scan point, 1806. Similarly to positioning the gates 
for the calculation and storing of the 2t data set file above, and 
prior to selecting the position and width of the gates for the 
front substrate echo (M M ) wavefonn and the back surface 
(BS) waveform echo, the user may canvas the fused file over 
the extent of the sample under evaluation to assist the user in 
deciding how wide (in time) the gates should be in time. 

Next in this example of the process, the step of using the 
stored values of 2 t and At in accordance with the algorithm to 
determine the velocity of the terahertz electromagnetic radia- 
tion in the sample on a scan point by scan point basis is 
performed 1807. The algorithm includes the step of determin- 
ing and storing, on a scan point by scan point basis, the 
velocity, V, of terahertz electromagnetic radiation in the 
sample without prior knowledge of the thickness. The veloc- 
ity, V, is determined by solving the equation, V=c(l-At/2x), 
for each the scan point and storing the determined velocity 
value in a velocity computer file on a scan point by scan point 
basis 1808. 

Similarly, using the stored values of 2x and At in the algo- 
rithm to determine the thickness of the sample on a scan point 
by scan point basis in accordance with the algorithm to deter- 
mine on a scan point by scan point basis, the thickness, d, of 
the sample without prior knowledge of the velocity of the 
terahertz electromagnetic radiation in the sample. The thick- 
ness is determined by solving the equation, d=c(2x-At)/2, for 
each the scan point and storing the thickness value in a thick- 
ness computer file 1809. 

Another exemplary process for non-destructive evaluation 
of a sample using a computer, a display and a computer 
program implementing an algorithm is disclosed in a block 
diagram 1700 in FIGS. 17-17C. The algorithm includes deter- 
mining the velocity, V, of terahertz electromagnetic radiation 
in the sample without prior knowledge of the thickness of the 
sample and determining the thickness of the sample without 
prior knowledge of the velocity of the terahertz electromag- 
netic radiation in the sample. The terahertz electromagnetic 
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radiation produced by a source spaced apart from the sample 
under test and propagated at the speed of light, c, in a medium 
located between the source and the sample. 

The exemplary process of FIGS. 17-17C includes the steps 
5 of: emitting and scanning terahertz electromagnetic radiation 
from the source, the terahertz radiation applied to a plurality 
of scan points arranged according to x and y coordinates of 
the sample, the plurality of scan points includes an area at 
least as large as the surface of the sample, 1701. Recording 
to and storing in a substrate echo waveform data set file, on a 
scan point by scan point basis, the substrate echo (M M ) wave- 
form of the terahertz electromagnetic radiation received from 
the substrate without the sample present is performed, 1702. 
Next, placing the sample on the substrate, 1703 is required so 
15 as to obtain front surface (FS) echo waveforms and back 
surface (BS) echo waveforms. Recording and storing in a 
back surface echo waveform data set file, on a scan point by 
scan point basis, the back surface echo (BS) waveform of the 
terahertz electromagnetic radiation received from the sub- 
20 strate with the sample present, 1704 is the next step per- 
formed. Further, the step of recording and storing in a front 
surface echo waveform data set file, on a scan point by scan 
point basis, the front surface echo (FS) waveform of the 
terahertz electromagnetic radiation received from the front 
25 surface (FS) of the sample, 1705 is performed. 

In the exemplary process of FIGS. 17-17C, the next step is 
recalling the stored substrate echo (M") waveform data set 
file, recalling the stored back surface echo (BS) waveform 
data set file, and recalling the stored front surface echo (FS) 
30 waveform data set file, 1706. Next, the step of opening and 
viewing, interactively and selectively, the substrate echo (M") 
waveform, the front surface echo (FS) waveform, and the 
back surface echo (BS) waveform, for a single scan point, the 
substrate echo (M M ) waveform, the front surface echo (FS) 
35 waveform, and the back surface echo (BS) waveform being 
superimposed in a display indicating amplitude and time base 
of the waveform, 1707 is performed. 

Next, viewing, interactively and selectively, an algorithm 
enabled display of all scan points and associated substrate 
40 (M M ) waveforms, associated front surface (FS) waveforms 
and associated back surface echo (BS) waveforms, for 
selected scan points of the sample arranged according to x and 
y scan coordinates, 1708 is performed. Gating, interactively, 
the front surface echo (FS) waveform and the substrate echo 
45 (M M ) waveform based on ranges determined from viewing 
selected scan points of the sample. Gating filters the wave- 
form portions beyond the time width of the gates, 1709. Next, 
the step of conditioning, the front surface echo (FS) wave- 
form, 1710 is performed. Optional steps of conditioning the 
50 front surface (FS) echo waveform may include processes 
selected from the group of wavelet denoising, wavelet recon- 
struction, and filtering, 1710 A. 

Next, the step of fusing the stored substrate echo (M M ) 
waveform data set file and the stored front surface echo (FS) 
55 waveform data set file into the stored back surface echo (BS) 
waveform data set file and producing a fused data file for each 
of the scan points, 1711 is performed. Once fused, the step of 
viewing and displaying the fused substrate (M M ) waveforms, 
front surface (FS) waveforms and associated back surface 
60 echo (BS) waveforms, interactively and selectively, for 
selected fused waveforms having x and y scan coordinates, 
according to an algorithm enabled display of all scan points, 
1712 is performed. This interactively allows the user to view 
the fused file set according to a specific data point and that 
65 data point coordinate is viewable by the user. 

Next, the exemplary process of FIGS. 17-17C includes the 
user selecting one of the fused data files, and interactively 



US 7,933,027 B1 


27 

gating the front surface (FS) echo waveform and interactively 
gating the back surface (BS) echo waveform based on ranges 
determined from viewing the selected fused waveforms, 

1713 . The user selects the gates as aforestated to calculate the 
difference in time, 2x, between the front surface (FS) echo 
waveforms and the back surface (BS) echo waveforms on a 
scan point by scan point basis generating and storing a 2r file, 

1714. The step of calculating the difference in time, 2x, 
between the front surface (FS) echo waveforms and the back 
surface (BS) echo waveforms on a scan point by scan point 
basis generating and storing a 2 t file is performed according 
to the algorithm by cross -correlating the waveforms on a scan 
point by scan point basis 1714A. 

Still referring to the process of FIGS. 17-17C, the step of 
calculating according to the algorithm the difference in time, 
At, between the front surface (BS) echo waveforms and the 
substrate (M M ) echo waveforms on a scan point by scan point 
basis generating and storing a At file is optionally performed 
according to the algorithm by cross -correlating the wave- 
forms on a scan point by scan point basis, 1714B. 

Still referring to the process of FIGS. 17-17C, the step of 
interactively gating the back surface (BS) echo waveform and 
interactively gating the substrate (M M ) echo waveform of the 
selected fused waveform file based on ranges determined 
from viewing the selected fused waveforms, 1715 is per- 
formed. Next, the step of calculating according to the algo- 
rithm the difference in time, At, between the front surface 
(BS) echo waveforms and the substrate (M") echo waveforms 
on a scan point by scan point basis generating and storing a At 
file, 1716 is performed. The step of calculating the difference 
in time, 2 t, between the front surface (FS) echo waveforms 
and the back surface (BS) echo waveforms on a scan point by 
scan point basis generating and storing a 2 t file is performed 
by comparing the peaks of the waveforms on a scan point by 
scan point basis, 1716A. Still referring to the exemplary 
process of FIGS. 17-17C, the step of calculating according to 
the algorithm the difference in time, At, between the front 
surface (BS) echo wavefonns and the substrate (M M ) echo 
waveforms on a scan point by scan point basis generating and 
storing a At file is optionally performed by comparing the 
peaks of the waveforms on a scan point by scan point basis, 
1716B. 

Still referring to the exemplary process of FIGS. 17-17C, 
the next step is determining according to the algorithm and 
storing, on a scan point by scan point basis, the velocity of 
terahertz electromagnetic radiation in the sample without 
prior knowledge of the thickness of the sample 1717. Next the 
results are displayed, on a scan point by scan point basis, the 
velocity of terahertz electromagnetic radiation in the sample 
without prior knowledge of the thickness of the sample, 
1717A. If the results are preferred in terms of density, the step 
of calculating, on a scan point by scan point basis, the density 
of the sample without prior knowledge of the thickness of the 
sample, 1717B is performed according to a look-up table 
based on known properties of the test sample under evalua- 
tion. 

According to the exemplary process of FIGS. 17-17C, the 
next step of determining and storing according to an algo- 
rithm, on a scan point by scan point basis, the velocity of 
terahertz electromagnetic radiation in the sample without 
prior knowledge of the thickness, is performed according to 
the algorithm solving for V using the equation: V=c(l -At/2x), 
1717C. Next, the step of determining according to the algo- 
rithm and storing, on a scan point by scan point basis, the 
thickness of the sample without prior knowledge of the veloc- 
ity of the terahertz electromagnetic radiation in the sample, 
1718 is performed. Once the velocity is determined it is 
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displayed, on a scan point by scan point basis, the thickness of 
the sample without prior knowledge of the velocity of the 
terahertz electromagnetic radiation in the sample, 1718A. 

And finally, according to the exemplary process set forth in 
5 FIGS. 17-17C, the step of determining according to the algo- 
rithm and storing, on a scan point by scan point basis, the 
thickness of the sample without prior knowledge of the veloc- 
ity of the terahertz electromagnetic radiation in the sample, is 
performed according to the algorithm solving for “d” using 
to the equation: d=c(2x-At)/2, 1718C. 

The step of fusing in the exemplary processes of FIGS. 
17-17C and 18-18A includes doubling the time base of the 
fused data waveform, 1790. 

The steps of gating and conditioning the substrate echo 
15 (M M ) and the front surface (FS) echo waveforms in the exem- 
plary process of FIGS. 17 and 18 include the step of delaying 
the substrate echo in time to better gate the substrate echo and 
then remove the delay when the step of fusing the substrate 
echo (M M ) waveform and the front side echo (FS) waveform 
20 into the back side echo waveform (BS) is performed, 1791. 

Another exemplary process for non-destructive evaluation 
of a sample using a computer, a display and a computer 
program implementing an algorithm is disclosed in a block 
diagram 1900 in FIG. 19. The process of FIG. 19 includes 
25 emitting and scanning terahertz electromagnetic radiation 
from said source, said terahertz radiation applied to a plurality 
of scan points arranged according to x and y coordinates of 
said source, said plurality of scan points includes an area at 
least as large as said surface of said sample, 1901. The process 
30 includes forming substrate (M M ) echo waveform files for each 
scan point 1902 and forming front surface (FS) echo wave- 
form files and back surface (BS) echo waveform files for each 
scan point 1903. 

According to the process of the example of FIG. 19, the 
35 next step includes superimposing, using a graphical user 
interface, said front surface (FS) echo waveform files, said 
back surface (BS) echo waveform files and said substrate 
echo (M M ) waveform for each scan point 1904. Next, the step 
of outputting and displaying an amplitude based image of one 
40 of said back surface (BS) echo files for inspection of said 
superimposed files 1905 is performed. 

Interactively, through a graphical user interface, inspecting 
said superimposed files; and selecting one of said superim- 
posed files for gating is performed 1906. Gating, using a 
45 graphical user interface, said front surface (FS) echo wave- 
form and said substrate echo (M") waveform and fusing said 
front surface (FS) echo waveform and said substrate echo 
(M M ) waveform into said back surface (BS) echo waveform 
on a scan point by scan point basis and outputting a composite 
50 image of said fused waveforms based on said amplitude of 
said back surface (BS) waveform 1907 is performed next. 
Selecting, using a graphical user interface, one of said fused 
waveforms and gating said waveform for determining 2x and 
At on a scan point by scan point basis, 1908 is performed. 
55 Next, using the algorithm, determining and storing, based on 
2x and At, on a scan point by scan point basis, the velocity of 
terahertz electromagnetic radiation in said sample without 
prior knowledge of the thickness of said sample is performed. 
Finally, the step of determining and storing, based on 2x and 
60 At, on a scan point by scan point basis, said thickness of said 
sample without prior knowledge of said velocity of said tera- 
hertz electromagnetic radiation in said sample 1909 is per- 
formed. 

FIG. 20 is an oscilloscope-like trace 2000 of a front surface 
65 (FS) echo waveform and a back surface (BS) echo waveform 
for the x=86, y=84 scan point coordinate. Reference numeral 
2001 denotes back surface (BS) echo waveform for the x=86, 
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y=84 scan point coordinate. Similarly, reference numeral 

2002 denotes the front surface (FS) echo waveform for the 
x=86, y-84 scan point coordinate. 

FIG. 20 A is an oscilloscope-like trace 2000A of a substrate 
echo (M M ) waveform for the x=86, y=84 scan point coordi- 
nate. Reference numeral 2003 denotes the substrate echo 
(M M ) waveform for the x=86, y=84 scan point coordinate. 

FIG. 20B is an oscilloscope-like trace 2000B of the front 
surface (FS) echo waveform 2002 and the substrate echo (M") 

2003 waveform superimposed into the back surface (BS) 
echo waveform 2001 for the x=86, y=84 scan point coordi- 
nate. 

FIG. 20C is an image 2000C of the back surface (BS) echo 
waveform illustrating the x-y scan points ranging from x=l 7 1 
to y=l 7 1 allowing the user to canvas any of the scan points to 
produce superimposed files having the front surface (FS) 
echo waveform, the substrate echo (M M ) waveform superim- 
posed into and onto the back surface (BS) echo waveform. 
Reference numeral 2004 denotes the x=86 scan coordinate, 
reference numeral 2005 denotes the y=84 scan coordinate, 
reference numeral 2006 denotes the cursor location/red circle 
and reference numeral 2007 indicates the gray scale ampli- 
tude meter. As such the present cursor position yields the 
information for the oscilloscope like traces illustrated in 
FIGS. 20, 20Aand20B. 

FIG. 20D is an image 2000D of the superimposed file for 
scan coordinates of x=20, y=l 50 and the corresponding oscil- 
loscope-like trace of the superimposed files having the front 
surface (FS) echo waveform, the substrate echo (M M ) wave- 
form superimposed into and onto the back surface (BS) echo 
waveform. Reference numeral 2008 illustrate the correspond- 
ing oscilloscope trace for cursor location x=20, y=150. 

FIG. 20E is an image 2000E of the superimposed file for 
scan coordinates of x=86, y=84 and the corresponding oscil- 
loscope-like trace of the superimposed files having the front 
surface (FS) echo waveform, the substrate echo (M M ) wave- 
form superimposed into and onto the back surface (BS) echo 
waveform. Reference numeral 2009 indicates the cursor loca- 
tion x=86, y=84 and the corresponding oscilloscope trace. 

FIG. 20F is an image 2000F of the superimposed file for 
scan coordinates of x=l 50, y=20 and the corresponding oscil- 
loscope-like trace of the superimposed files having the front 
surface (FS) echo waveform, the substrate echo (M M ) wave- 
form superimposed into and onto the back surface (BS) echo 
waveform. Reference numeral 2010 illustrates the cursor 
location x=150, y=20 and the corresponding oscilloscope- 
like trace. 

FIG. 20G is a screen shot 2000G from the software pro- 
gram implementing signal condition of the front surface (FS) 
echo waveform and delaying, by a specified amount of time, 
the substrate echo (M M ) waveform to facilitate gating of the 
substrate echo (M M ) waveform. Reference numeral 2011 
denotes a time shift input dialog box where the operator 
inserts a time delay to separate the substrate echo (M M ) from 
the back surface (BS) echo waveform. Reference numeral 
2003A indicates the position of the substrate echo which has 
been time shifted to separate it from the back surface (BS) 
echo waveform. 

FIG. 20H is an oscilloscope-like trace 2000H of the time 
shifted substrate echo (M M ) waveform and gating of the front 
surface (FS) echo waveform and the time shifted substrate 
(M M ) waveform. Reference numeral 2012 denotes the front 
surface (FS) echo waveform gate and reference numeral 2013 
indicates the substrate echo (M M ) waveform gate. These gates 
are placed as stated prior to fusing the waves as is illustrated 
in FIG. 201. 
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FIG. 201 is a screen shot 20001 similar to FIG. 20G from the 
software program implementing signal condition of the front 
surface (FS) echo waveform and delaying, by a specified 
amount of time, the substrate echo (M M ) waveform to facili- 
5 tate gating of the substrate echo (M M ) waveform and in addi- 
tion indicating the “fuse gated waves” boolean button. Ref- 
erence numeral 2014 globally denotes the boolean buttons for 
front surface (FS) echo wave form signal condition selectable 
options and reference numeral 2015 is the boolean button for 
to fusing the gated waves. 

FIG. 20J is an oscilloscope-like trace 2000J of the fused 
waveform 2016 containing the front surface (FS) echo, the 
substrate echo (M M ) and the back surface (BS) echo. FIG. 20K 
is the oscilloscope-like trace 2000K of the fused waveform 
15 containing the front surface (FS) echo, the substrate echo 
(M M ) and the back surface (BS) echo of FIG. 20J illustrating 
gates 2017, 2018 placed over the front surface (FS) echo and 
the back surface (BS) echo, respectively. 

FIG. 20L is the oscilloscope-like trace 2000L of the fused 
20 waveform containing the front surface (FS) echo, the sub- 
strate echo (M M ) and the back surface (BS) echo of FIG. 20 J 
illustrating gates 2019, 2020 placed over the substrate echo 
(M M ) and the back surface (BS) echo, respectively. 

FIG. 20M is a computer generated image 2000M of the 
25 thickness-independent velocity image generated by the com- 
puter program. Reference numeral 2021 is the velocity ampli- 
tude gray scale indicator associated with the image of FIG. 
20. Reference numeral 2022 denotes the cursor position of the 
velocity image. Check Box 2025 for converting velocity 
30 image to a density image is indicated in FIG. 20M enabling 
the user to convert from the velocity image to the density 
image of FIG. 20N. 

FIG. 2 ON is a computer generated image 2 00 ON of the 
density image corresponding to the velocity image of FIG. 
35 2 0M generated by the computer program. Reference numeral 
2023 indicates the density amplitude gray scale indicator and 
reference numeral 2024 indicates the cursor position of the 
density image. In both FIGS. 20M and 20N enable the user to 
canvas the images to obtain site specific information for the 
40 scan coordinates. 

FIG. 21 is a screen shot 2100 of the program illustrating a 
histogram of the density plot and a density image after the 
outliers have been removed. The remove outliers button is 
illustrated in the screen shot. The screen shot 2100 is obtained 
45 from the density image wherein processing is selected fol- 
lowed by “advanced” and remove outliers. Reference 
numeral denotes 2101 is the minimum value of acceptable 
density data points and reference numeral 2102 indicates the 
maximum value of acceptable density data points. Outliers 
50 2103 are indicated in the histogram and button 2104 prompts 
the program to remove outliers. Reference numeral 2105 
denotes the true density image. 

REFERENCE NUMERALS 

55 

100 — schematic view of the terahertz measuring system 

101 — transceiver 

102 — terahertz source 

103 — terahertz receiver 

60 1 04 — front surface of metal substrate 
104 A — metal substrate 

105 — back surface of metal substrate 

106 — beginning of void 

107 — end of void/front surface of metal substrate 
65 1 08 — foam, silicon nitride or other dielectric 

109 — front surface of foam, silicon nitride or other dielectric 

110 — gate for signal analysis 
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200 — density of sprayed on foam insulation versus velocity 
of terahertz electromagnetic radiation therein 

200A — density of silicon nitride versus velocity of terahertz 
electromagnetic radiation therein 

300 — schematic view of terahertz measuring system and 5 
graph of output voltages of respective signals versus time 

301 — thin dielectric sheet 

400 — graph of uncertainty (in percent) of thickness indepen- 
dent velocities as a function of 2 t, At and V. 

500 — graph of uncertainty (in percent) of velocity as a func- 
tion of thickness, d. 

600 — schematic of step wedge foam blocks 

601 — aluminum substrate 

602 — first foam block 15 

603 — second foam block 

604 — third foam block 

605 — fourth foam block 

700 — plot of 2 t as a function of thickness and density 

70 0A — schematic illustration of thickness variation 20 

70 0B — schematic illustration of density variation 

800 — physically measured density map for a 6 by 15 set of 
foam blocks 

800A — density plot 

801 — comparison portion of physically measured density 25 
map 

802 — comparison portion of physically measured density 
map 

900 — terahertz density map 

900A — terahertz density plot 30 

901 — comparison portion of thickness independent terahertz 
density map 

902 — comparison portion of thickness independent terahertz 
density map 

1000 — hand measured thickness map for a 6 by 1 5 set of foam 35 
blocks 

10 00 A — illustrates the thickness in centimeters by shade of 
gray 

1001 — ellipse indicating an area of the hand measured thick- 
ness map to be compared to a terahertz thickness map 40 

1100 — a terahertz thickness image for the 6 by 15 set of foam 
blocks (determined independently of velocity) according 
to the methodology of the invention 

11 00 A — thickness by shade of gray in centimeters for the 
terahertz thickness map of FIG. 11 45 

1101 — ellipse indicating an area of the terahertz thickness 
map to be compared to a physically measured thickness 
map 

1200 — schematic diagram of a process for measuring the 
velocity of terahertz electromagnetic radiation in a mate- 50 
rial sample without prior knowledge of the thickness of the 
sample, the terahertz electromagnetic radiation produced 
by a source spaced apart from the sample under test at the 
speed of light, c, in a medium located between the source 
and the sample, the sample residing on the substrate 55 

1201 — emitting terahertz electromagnetic radiation from the 
source 

1202 — measuring the travel time of the terahertz electromag- 

netic radiation to and from the substrate without the sample 
present, t" (t-double prime) 60 

1203 — placing the sample on the substrate 

1204 — measuring the travel time of the terahertz electromag- 
netic radiation to and from the substrate with the sample/ 
dielectric present, f (t-prime) 

1205 — subtracting t" from f to determine the transmission 65 
time difference with the sample present and without the 
sample present, At 
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1206 — measuring the travel time of the terahertz electromag- 
netic radiation to and from the sample, 2 t l 

1207 — determining the time difference between the FS echo 
and the BS echo, 2x 

1208 — placing a sheet of plastic paper onto the sample to 
create an adequate dielectric mismatch 

1209 — dividing At by 2x and determining the quotient 

1210 — subtracting the quotient from 1 to obtain a factor 

1211 — multiplying the factor to determine the velocity of the 
terahertz electromagnetic radiation in the sample 

1212 — determining the micro structural variation of the 
sample 

1213 — determining the density of the material 

1214 — evaluating, in a plurality of locations, the sample for 
micro structural variations; and mapping the variations by 
location 

1215 — controlling the spacing between the source and the 
substrate 

1216 — determining, initially, the approximate time location 
of the sample from the terahertz radiation source 

1300 — schematic diagram of a process for measuring the 
thickness of a material sample using terahertz electromag- 
netic radiation in a material sample without prior knowl- 
edge of the velocity of the terahertz electromagnetic radia- 
tion in the sample, the terahertz electromagnetic radiation 
produced by a source spaced apart from the sample under 
test and propagated at the speed of light, c, in a medium 
located between the source and the sample, the sample 
residing on the substrate 

1301 — emitting terahertz electromagnetic radiation from the 
source; 

1302 — measuring the travel time of the terahertz electromag- 
netic radiation to and from the substrate without the sample 
present, t” (t-double prime) 

1303 — placing the sample on the substrate 

1304 — measuring the travel time of the terahertz electromag- 
netic radiation to and from the substrate with the sample 
present, t' (t-prime) 

1305 — subtracting t" from f to determine the transmission 
time difference with the sample present and without the 
sample present, At 

1306 — measuring the travel time of the terahertz electromag- 
netic radiation to and from the sample, 2 t l 

1307 — determining the time difference between the FS echo 
and the BS echo, 2 t 

1308 — subtracting At from 2 t to obtain a subtraction result 

1309 — dividing the subtraction result by 2 to obtain a quo- 
tient 

1310 — multiplying the quotient by, c, to obtain the thickness 
of the sample 

1312 — evaluating, in a plurality of locations, the sample for 
thickness variations; and mapping the thickness variations 
by location 

1400 — schematic diagram of a process for measuring the 
thickness of a material sample using terahertz electromag- 
netic radiation in a material sample without prior knowl- 
edge of the velocity of the terahertz electromagnetic radia- 
tion in the sample, the terahertz electromagnetic radiation 
produced by a source spaced apart from the sample under 
test and propagated at the speed of light, c, in a medium 
located between the source and the sample under test, the 
sample residing on the substrate 

1401 — emitting terahertz electromagnetic radiation from the 
source 

1402 — measuring the travel time of the terahertz electromag- 
netic radiation to and from the substrate with the sample 
present, t' (t-prime) 
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1403 — placing the sample on the substrate 

1404 — measuring the travel time of the terahertz electromag- 
netic radiation to and from the substrate without the sample 
present, t" (t-double prime) 

1405 — subtracting t" from f to determine the transmission 
time difference with the sample present and without the 
sample present, At 

1406 — measuring the travel time of the terahertz electromag- 
netic radiation to and from the sample, 2 1 1 

1407 — determining the time difference between the FS echo 
and the BS echo, 2x 1408 — dividing At by 2 t and deter- 
mining the quotient 

1409 — subtracting the quotient from 1 to obtain a factor 

1410 — multiplying the factor to determine the velocity of the 
terahertz electromagnetic radiation in the sample 

1412 — subtracting At from 2x to obtain a subtraction result 

1413 — dividing the subtraction result by 2 to obtain a quo- 
tient 

1414 — multiplying the quotient by, c, to obtain the thickness 
of the sample 

1416 — evaluating, in a plurality of locations, the sample for 
micro structural variations and for thickness variations, and 
mapping the micro structural and thickness variations by 
location. 

1500 — fused waveform from a “fused” data file produced by 
merging the FS, BS, and M" scan data sets for a set of foam 
blocks 

1600 — another process for simultaneously measuring the 
velocity of terahertz electromagnetic radiation in a mate- 
rial sample without prior knowledge of the thickness of the 
sample and for measuring the thickness of a material 
sample using terahertz electromagnetic radiation in a mate- 
rial sample without prior knowledge of the velocity of the 
terahertz electromagnetic radiation in the sample, the tera- 
hertz electromagnetic radiation produced by a source 
spaced apart from the substrate and propagated at the speed 
of light, c, in a medium located between the source and the 
substrate. 

1601 — merging the FS, BS, and M" scan data sets into a fused 
data file. 

1602 — producing a “fused” waveform 1602. 

1603 — calculating time delays between echos. 

1604 — determining precision thickness-independent veloc- 
ity images that map microstructure. 

1605 — determining precision micro structure-independent 
thickness images that map thickness. 

1606 — calculating the time delays and determining precision 
thickness and microstructure are performed by the step of 
cross-correlating entire waveforms FS, BS and M M . 

1607 — calculating the time delays and determining precision 
thickness and microstructure are performed by precisely 
identifying peaks of FS, BS and M" signals. 

1608 — gating and processing the FS signal prior to merging 
the data into a fused data file to account for low FS signal 
(if needed) 

1609 — amplifying the processed FS signal prior to merging 
the data into a fused data file to account for FS signal (if 
needed) 

1610 — applying a dielectric sheet on the front surface, fol- 
lowed by processing and gating the FS signal prior to 
merging the data into a fused data file to account for low FS 
signal (if needed) 

1700 — schematic of exemplary computer implemented pro- 
cess 

1701 — emitting and scanning terahertz electromagnetic 
radiation from the source, the terahertz radiation applied to 
a plurality of scan points arranged according to x and y 
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coordinates of the sample, the plurality of scan points 
includes an area at least as large as the surface of the sample 

1702 — recording and storing in a substrate echo waveform 
data set file, on a scan point by scan point basis, the sub- 

5 strate echo (M") waveform of the terahertz electromag- 
netic radiation received from the substrate without the 
sample present; 

1703 — placing the sample on the substrate 

1704 — recording and storing in a back surface echo wave- 

10 form data set file, on a scan point by scan point basis, the 

back surface echo (BS) waveform of the terahertz electro- 
magnetic radiation received from the substrate with the 
sample present 

15 1705 — recording and storing in a front surface echo wave- 
form data set file, on a scan point by scan point basis, the 
front surface echo (FS) waveform of the terahertz electro- 
magnetic radiation received from the front surface (FS) of 
the sample; 

20 1 706 — recalling the stored substrate echo (M M ) waveform 
data set file, recalling the stored back surface echo (BS) 
waveform data set file, and recalling the stored front sur- 
face echo (FS) waveform data set file 

1707 — opening and viewing, interactively and selectively, 

25 the substrate echo (M") waveform, the front surface echo 

(FS) waveform, and the back surface echo (BS) waveform, 
for a single scan point, the substrate echo (M M ) waveform, 
the front surface echo (FS) waveform, and the back surface 
echo (BS) waveform being superimposed in a display indi- 

30 eating amplitude and time base of the waveform 

1708 — viewing, interactively and selectively, an algorithm 
generated display of all scan points and associated sub- 
strate (M M ) waveforms, associated front surface (FS) wave- 
forms and associated back surface echo (BS) waveforms, 

35 for selected scan points of the sample arranged according 

to x and y scan coordinates 

1709 — gating, interactively, the front surface echo (FS) 
waveform and the substrate echo (M") waveform based on 
ranges determined from viewing selected scan points of the 

40 sample, the gating filters the scan points having waveforms 
beyond the time width of the gates; 

1710 — conditioning, the front surface echo (FS) waveform 

1710 A — optionally conditioning the front surface (FS) echo 

waveform includes processes selected from the group of 

45 wavelet denoising, wavelet reconstruction, and filtering 

1711 — fusing the stored substrate echo (M M ) waveform data 
set file and the stored front surface echo (FS) waveform 
data set file into the stored back surface echo (BS) wave- 
form data set file and producing a fused data file for each of 

50 the scan points; 

1712 — viewing and displaying the fused substrate (M M ) 
waveforms, front surface (FS) waveforms and associated 
back surface echo (BS) waveforms, interactively and selec- 
tively, selected x and y scan coordinates, according to an 

55 algorithm generated display of all scan points; 

1713 — selecting one of the fused data files, and interactively 
gating the front surface (FS) echo waveform and interac- 
tively gating the back surface (BS) echo waveform based 
on ranges determined from viewing the selected fused 

60 waveforms; 

1714 — calculating the difference in time, 2x, between the 
front surface (FS) echo waveforms and the back surface 
(BS) echo waveforms on a scan point by scan point basis 
generating and storing a 2x file; 

65 171 4 A — calculating the difference in time, 2x, between the 
front surface (FS) echo waveforms and the back surface 
(BS) echo waveforms on a scan point by scan point basis 
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generating and storing a 2x file is performed according to 
the algorithm by cross-correlating the waveforms on a scan 
point by scan point basis 

1714B — optionally calculating the difference in time, 2 t, 
between the front surface (FS) echo waveforms and the 
back surface (BS) echo waveforms on a scan point by scan 
point basis generating and storing a 2x file is performed 
according the algorithm by comparing the peaks of the 
waveforms on a scan point by scan point basis 

1715 — interactively gating the back surface (BS) echo wave- 
form and interactively gating the substrate (M M ) echo 
waveform of the selected fused waveform file based on 
ranges determined from viewing the selected fused wave- 
forms; 

1716 — calculating according to the algorithm the difference 
in time, At, between the front surface (BS) echo waveforms 
and the substrate (M M ) echo waveforms on a scan point by 
scan point basis generating and storing a At file 

171 6 A — optionally calculating according to the algorithm 
the difference in time, At, between the front surface (BS) 
echo waveforms and the substrate (M M ) echo waveforms on 
a scan point by scan point basis generating and storing a At 
file is performed according to the algorithm by cross-cor- 
relating the waveforms on a scan point by scan point basis 

1716B — optionally calculating according to the algorithm 
the difference in time, At, between the front surface (BS) 
echo waveforms and the substrate (M") echo waveforms on 
a scan point by scan point basis generating and storing a At 
file is performed according to the algorithm by comparing 
the peaks of the waveforms on a scan point by scan point 
basis 

1717 — determining according to the algorithm and storing, 
on a scan point by scan point basis, the velocity of terahertz 
electromagnetic radiation in the sample without prior 
knowledge of the thickness of the sample 

1717A — displaying, on a scan point by scan point basis, the 
velocity of terahertz electromagnetic radiation in the 
sample without prior knowledge of the thickness of the 
sample 

1717B — calculating, on a scan point by scan point basis, the 
density of in the sample without prior knowledge of the 
thickness of the sample 

1717C — determining and storing according to an algorithm, 
on a scan point by scan point basis, the velocity of terahertz 
electromagnetic radiation in the sample without prior 
knowledge of the thickness, according to the algorithm 
implementing the equation: V=c(l-At/2x) 

1718 — determining according to the algorithm and storing, 
on a scan point by scan point basis, the thickness of the 
sample without prior knowledge of the velocity of the 
terahertz electromagnetic radiation in the sample 

1718 A — displaying, on a scan point by scan point basis, the 
thickness of the sample without prior knowledge of the 
velocity of the terahertz electromagnetic radiation in the 
sample 

1718B — step of determining according to the algorithm and 
storing, on a scan point by scan point basis, the thickness of 
the sample without prior knowledge of the velocity of the 
terahertz electromagnetic radiation in the sample, is per- 
formed according to the algorithm implementing the equa- 
tion: d=c(2x-At)/2 

1790 — step of fusing in the exemplary processes of FIGS. 17 
and 18 includes doubling the time base of the fused data 
waveform. 

1791 — steps of gating and conditioning the substrate echo 
(M M ) and the front surface (FS) echo waveforms in the 
exemplary process of FIGS. 17 and 18 include the step of 
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delaying the substrate echo in time to better gate the sub- 
strate echo and then remove the delay when the step of 
fusing the substrate echo (M M ) waveform and the front side 
echo (FS) waveform into the back side echo waveform 
5 (BS) is performed. 

1800 — schematic of exemplary computer implemented pro- 
cess 

1801 — emitting and scanning terahertz electromagnetic 
radiation from the source, the terahertz radiation applied to 

to a plurality of scan points arranged according to x and y 
coordinates of the scanning, the plurality of scan points 
includes an area at least as large as the surface of the 
sample, 1801. 

1802 — recording and storing in a substrate echo waveform 
15 data set file, on a scan point by scan point basis, a substrate 

echo (M M ) waveform of the terahertz electromagnetic 
radiation received from the substrate without the sample 
present, and recording and storing front surface (FS) wave- 
form and back surface (BS) waveform of the terahertz 
20 electromagnetic radiation received with the sample present 

1803 — superimposing the substrate echo (M M ) waveform, the 
front surface (FS) echo waveform, and the back surface 
(BS) echo waveform together 

1804 — fusing the substrate echo (M M ) waveform and the front 
25 surface echo (FS) waveform into the back surface (BS) 

echo waveform 

1805 — gating the echo waveforms and determine the time 
difference between the front surface echo (FS) waveform 
and the back surface (BS) echo waveform, 2x, and storing 

30 values of 2x for each scan point 

1806 — interactively gating the echo waveforms and deter- 
mining the time difference between the substrate echo 
wave (M M ) waveform and the back surface (BS) waveform, 
At, and storing values of At for each scan point 

35 1807 — using the stored values of 2x and At in accordance 
with the algorithm to determine the velocity of the terahertz 
electromagnetic radiation in the sample on a scan point by 
scan point basis is performed 

1808 — determining and storing, on a scan point by scan point 

40 basis, the velocity, V, of terahertz electromagnetic radiation 

in the sample without prior knowledge of the thickness, 
velocity, V, is determined by solving the equation, V=c(l - 
At/2x), for each the scan point and storing the determined 
velocity value in a velocity computer file on a scan point by 
45 scan point basis 

1809 — using the stored values of 2x and At in the algorithm to 
determine the thickness of the sample on a scan point by 
scan point basis in accordance with the algorithm to deter- 
mine on a scan point by scan point basis, the thickness, d, 

50 of the sample without prior knowledge of the velocity of 
the terahertz electromagnetic radiation in the sample, 
thickness is determined by solving the equation, d=c(2x- 
At)/2, for each the scan point and storing the thickness 
value in a thickness computer file 
55 1900 — schematic of exemplary computer implemented pro- 
cess 

1901 — emitting and scanning terahertz electromagnetic 
radiation from said source, the terahertz radiation applied 
to a plurality of scan points arranged according to x and y 

60 coordinates of the source, the plurality of scan points 
includes an area at least as lai^e as the surface of the sample 

1902 — forming a substrate (M H ) echo waveform files for each 
scan point 

1903 — forming a front surface (FS) echo waveform file and a 
65 back surface (BS) echo waveform file for each scan point 

1904 — superimposing, using a graphical user interface, said 
front surface (FS) echo waveform files, the back surface 
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(BS) echo waveform files and the substrate echo (M") 
waveform for each scan point 

1905 — outputting and displaying an amplitude based image 
of one of the superimposed files for inspection of the super- 
imposed files 

1906 — inspecting said superimposed files and selecting one 
of said superimposed files for gating 

1907 — gating, using a graphical user interface, the front sur- 
face (FS) echo waveform and the substrate echo (M M ) 
waveform and fusing said front surface (FS) echo wave- 
form and the substrate echo (M M ) waveform into the back 
surface (BS) echo waveform on a scan point by scan point 
basis and outputting an image of said fused waveforms 

1908 — selecting, using a graphical user interface, one of said 
fused waveforms and gating said waveform for determin- 
ing 2x and At on a scan point by scan point basis 

1909 — determining and storing, based on 2x and At, on a scan 
point by scan point basis, the velocity of terahertz electro- 
magnetic radiation in the sample without prior knowledge 
of the thickness of the sample is performed and determin- 
ing and storing, based on 2x and At, on a scan point by scan 
point basis, the thickness of the sample without prior 
knowledge of the velocity of the terahertz electromagnetic 
radiation in the sample 

2000 — oscilloscope-like trace of a front surface (FS) echo 
waveform and a back surface (BS) echo waveform for the 
x=86, y-84 scan point coordinate 

2000A — oscilloscope-like trace of a substrate echo (M") 
waveform for the x=86, y-84 scan point coordinate 

2000B — oscilloscope-like trace of the front surface (FS) echo 
waveform and the substrate echo (M") waveform superim- 
posed into the back surface (BS) echo waveform for the 
x=86, y-84 scan point coordinate 

2000C — image of the back surface (BS) echo waveform illus- 
trating the x-y scan points ranging from x=171 to y=171 
allowing the user to canvas any of the scan points to pro- 
duce superimposed files having the front surface (FS) echo 
waveform, the substrate echo (M") waveform superim- 
posed into and onto the back surface (BS) echo waveform 

2000D — image of the superimposed file for scan coordinates 
of x=20, y=150 and the corresponding oscilloscope-like 
trace of the superimposed files having the front surface 
(FS) echo waveform, the substrate echo (M M ) waveform 
superimposed into and onto the back surface (BS) echo 
waveform 

2000E — image of the superimposed file for scan coordinates 
of x=86, y=84 and the corresponding oscilloscope-like 
trace of the superimposed files having the front surface 
(FS) echo waveform, the substrate echo (M") waveform 
superimposed into and onto the back surface (BS) echo 
waveform 

2000F — image of the superimposed file for scan coordinates 
of x=l 50, y=20 and the corresponding oscilloscope-like 
trace of the superimposed files having the front surface 
(FS) echo waveform, the substrate echo (M") waveform 
superimposed into and onto the back surface (BS) echo 
waveform 

2000G — screen shot from the software program implement- 
ing signal condition of the front surface (FS) echo wave- 
form and delaying, by a specified amount of time, the 
substrate echo (M H ) waveform to facilitate gating of the 
substrate echo (M M ) waveform 

2000F — oscilloscope-like trace of the time shifted substrate 
echo (M M ) waveform and gating of the front surface (FS) 
echo waveform and the time shifted substrate (M M ) wave- 
form 
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20001 — screen shot similar to FIG. 20G from the software 
program implementing signal condition of the front sur- 
face (FS) echo waveform and delaying, by a specified 
amount of time, the substrate echo (M") waveform to facili- 

5 tate gating of the substrate echo (M M ) waveform and in 
addition indicating the “fuse gated waves” boolean button 

2000J — oscilloscope-like trace of the fused waveform con- 
taining the front surface (FS) echo, the substrate echo (M M ) 
and the back surface (BS) echo 

to 2000K — oscilloscope-like trace of the fused waveform con- 
taining the front surface (FS) echo, the substrate echo (M M ) 
and the back surface (BS) echo of FIG. 20J illustrating 
gates placed over the front surface (FS) echo and the back 
surface (BS) echo 

15 2000L — oscilloscope-like trace of the fused waveform con- 
taining the front surface (FS) echo, the substrate echo (M M ) 
and the back surface (BS) echo of FIG. 20J illustrating 
gates placed over the substrate echo (M M ) and the back 
surface (BS) echo. 

20 2000M — computer generated image of the thickness -inde- 
pendent velocity image generated by the computer pro- 
gram. 

2 00 ON — computer generated image of the density image cor- 
responding to the velocity image of FIG. 2000M. 

25 2 001 — back surface (BS) echo waveform for the x=86, y=84 
scan point coordinate 

2002 — front surface (FS) echo waveform for the x=86, y=84 
scan point coordinate 

2003 — substrate echo (M") waveform for the x=86, y=84 

30 scan point coordinate 

2003 — time shifted substrate echo (M") waveform 

2004 — x=86 scan coordinate 

2005 — y=84 scan coordinate 

2006 — cursor location/red circle 

35 2007 — gray scale amplitude meter 

2008 — cursor location x=20, y=150 

2009 — cursor location x=86, y=84 

2010 — cursor location x=150, y=20 

2011 — time shift input dialog box 

40 2012 — front surface (FS) echo waveform gate 

2013 — substrate echo (M") waveform gate 

2014 — boolean buttons for front surface (FS) echo wave form 
signal condition selectable options 

2015 — boolean button for fusing gated waves 

45 2016 — fused waveform comprising front surface (FS) echo 
waveforms, substrate echo (M M ) waveforms, and back sur- 
face (BS) echo waveforms 

2017 — gate over the front surface (FS) echo waveform of the 
fused wavetrain 

50 2018 — gate over the back surface (BS) echo waveform of the 
fused wavetrain 

2019 — gate over the substrate echo (M H ) echo waveform of 
the fused wavetrain 

2020 — gate over the back surface (BS) echo waveform of the 

55 fused wavetrain 

2021 — velocity amplitude gray scale indicator 

2022 — cursor position of the velocity image 

2023 — density amplitude gray scale indicator 

2024 — cursor position of the density image 

60 2025 — check box for converting velocity image to a density 
image 

2100 — is a screen shot of the program illustrating a histogram 
of the density plot and a density image after the outliers 
have been removed. 

65 2101 — minimum value of acceptable density data points 

2102 — maximum value of acceptable density data points 

2103 — outliers 
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2104 — remove outlier button 

2105 — true density image 

BS — pulse that travels from the transceiver to the reflector 
plate and back to the transceiver with the sample present 
FS — pulse that travels from the transceiver to the front sur- 5 
face of the sample and back to the transceiver 
L — distance between transceiver and sample 
M M — pulse that travels from the transceiver to the reflector 
plate and back to the transceiver without the sample present 
c — speed of light in a medium located between the source of to 
the terahertz radiation and the substrate 
d — sample thickness 

f (t-prime) — travel time of the terahertz electromagnetic 
radiation to and from the substrate with the sample present; 
t" (t-double prime) — travel time of the terahertz electromag- 15 
netic radiation to and from the substrate without the sample 
present; 

At — transmission time difference (t' minus t") with the 
sample present and without the sample present; 

2t l — measuring the travel time of the terahertz electromag- 20 
netic radiation to and from the sample; 

2t — time difference between the FS echo and the BS echo, 2% 

V — velocity 

Those skilled in the art will readily recognize that the 
invention has been set forth by way of example only. Accord- 25 
ingly, those skilled in the art will recognize that changes may 
be made to the invention without departing from the spirit and 
scope of the attached claims. 

The invention claimed is: 

1. A process for non-destructive evaluation of a sample 30 
using a computer, a display and a computer program imple- 
menting an algorithm, said algorithm includes determining 
the velocity of terahertz electromagnetic radiation in said 
sample without prior knowledge of the thickness of said 
sample and determining said thickness of said sample without 35 
prior knowledge of said velocity of said terahertz electromag- 
netic radiation in said sample, said terahertz electromagnetic 
radiation produced by a source, namely, a transceiver, spaced 
apart from said sample under test and propagated at the speed 
of light, c, in a medium located between said source, namely, 40 
a transceiver, and said sample, comprising the steps of: 
emitting and scanning terahertz electromagnetic radiation 
from said source, namely, said transceiver, said terahertz 
radiation applied to a plurality of scan points arranged 
according to x and y coordinates of said sample, said 45 
plurality of scan points includes an area at least as large 
as said surface of said sample; 
recording and storing in a substrate echo waveform data set 
file, on a scan point by scan point basis, said substrate 
echo (M") waveform of said terahertz electromagnetic 50 
radiation received from said substrate without said 
sample present; 

placing said sample on said substrate; 
recording and storing in a back surface echo waveform data 
set file, on a scan point by scan point basis, the back 55 
surface echo (BS) waveform of said terahertz electro- 
magnetic radiation received from said substrate with 
said sample present; 

recording and storing in a front surface echo waveform 
data set file, on a scan point by scan point basis, the front 60 
surface echo (FS) waveform of said terahertz electro- 
magnetic radiation received from said front surface (FS) 
of the sample; 

recalling said stored substrate echo (M") waveform data set 
file, recalling said stored back surface echo (BS) wave- 65 
form data set file, and recalling said stored front surface 
echo (FS) waveform data set file; 
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opening and viewing, interactively and selectively, said 
substrate echo (M") waveform, said front surface echo 
(FS) waveform, and said back surface echo (BS) wave- 
form, for a single scan point, said substrate echo (M") 
waveform, said front surface echo (FS) waveform, and 
said back surface echo (BS) waveform being superim- 
posed in a display indicating amplitude and time base of 
said waveform; 

viewing, interactively and selectively, an algorithm gener- 
ated display of all scan points and associated substrate 
(M") waveforms, associated front surface (FS) wave- 
forms and associated back surface echo (BS) wave- 
forms, for selected scan points of said sample arranged 
according to x and y scan coordinates; 
gating, interactively, said front surface echo (FS) wave- 
form and said substrate echo (M") waveform based on 
ranges determined from viewing selected scan points of 
said sample, said gating filters said scan points having 
waveforms beyond the time width of said gates; 
conditioning, said front surface echo (FS) waveform; 
fusing said stored substrate echo (M M ) waveform data set 
file and said stored front surface echo (FS) waveform 
data set file into said stored back surface echo (BS) 
waveform data set file and producing a fused data file for 
each of said scan points; 

viewing and displaying said fused substrate (M") wave- 
forms, front surface (FS) waveforms and associated 
back surface echo (BS) waveforms, interactively and 
selectively, for selected x and y scan coordinates, 
according to an algorithm generated display of all scan 
points; 

selecting one of said fused data files, and interactively 
gating said front surface (FS) echo waveform and inter- 
actively gating said back surface (BS) echo waveform of 
one said selected fused waveform file based on ranges 
determined from viewing said previously fused wave- 
forms; 

calculating the difference in time, 2 t, between said front 
surface (FS) echo waveforms and said back surface (BS) 
echo waveforms on a scan point by scan point basis 
generating and storing a 2 t file; 
interactively gating said back surface (BS) echo waveform 
and interactively gating said substrate (M") echo wave- 
form of said selected fused waveform file based on 
ranges determined from viewing said selected fused 
waveforms; 

calculating according to said algorithm the difference in 
time, At, between said front surface (BS) echo wave- 
forms and said substrate (M H ) echo waveforms on a scan 
point by scan point basis generating and storing a At file; 
determining according to said algorithm and storing, on a 
scan point by scan point basis, the velocity of terahertz 
electromagnetic radiation in said sample without prior 
knowledge of the thickness of said sample; and, 
determining according to said algorithm and storing, on a 
scan point by scan point basis, said thickness of said 
sample without prior knowledge of said velocity of said 
terahertz electromagnetic radiation in said sample. 

2. A process for non-destructive evaluation of a sample 
using a computer, a display and a computer program imple- 
menting an algorithm, said algorithm includes determining 
the velocity of terahertz electromagnetic radiation in said 
sample without prior knowledge of the thickness of said 
sample and determining said thickness of said sample without 
prior knowledge of said velocity of said terahertz electromag- 
netic radiation in said sample, said terahertz electromagnetic 
radiation produced by a source, namely, a transceiver, spaced 
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apart from said sample under test and propagated at the speed 
of light, c, in a medium located between said source, namely, 
said transceiver, and said sample as claimed in claim 1, fur- 
ther comprising the steps of: 

displaying, on a scan point by scan point basis, said veloc- 
ity of terahertz electromagnetic radiation in said sample 
without prior knowledge of the thickness of said sample; 
and, 

displaying, on a scan point by scan point basis, said thick- 
ness of said sample without prior knowledge of said 
velocity of said terahertz electromagnetic radiation in 
said sample. 

3. A process for non-destructive evaluation of a sample 
using a computer, a display and a computer program imple- 
menting an algorithm, said algorithm includes determining 
the velocity of terahertz electromagnetic radiation in said 
sample without prior knowledge of the thickness of said 
sample and determining said thickness of said sample without 
prior knowledge of said velocity of said terahertz electromag- 
netic radiation in said sample, said terahertz electromagnetic 
radiation produced by a source, namely, a transceiver, spaced 
apart from said sample under test and propagated at the speed 
of light, c, in a medium located between said source, namely, 
said transceiver, and said sample as claimed in claim 1, fur- 
ther comprising the steps of: 

calculating, on a scan point by scan point basis, the density 
of terahertz electromagnetic radiation in said sample 
without prior knowledge of the thickness of said sample. 

4. A process for non-destructive evaluation of a sample 
using a computer, a display and a computer program imple- 
menting an algorithm, said algorithm includes determining 
the velocity of terahertz electromagnetic radiation in said 
sample without prior knowledge of the thickness of said 
sample and determining said thickness of said sample without 
prior knowledge of said velocity of said terahertz electromag- 
netic radiation in said sample, said terahertz electromagnetic 
radiation produced by a source, namely, a transceiver, spaced 
apart from said sample under test and propagated at the speed 
of light, c, in a medium located between said source, namely, 
said transceiver, and said sample as claimed in claim 3, fur- 
ther comprising the steps of: 

displaying, on a scan point by scan point basis, the density 
of terahertz electromagnetic radiation in said sample 
without prior knowledge of the thickness of said sample. 

5. A process for non-destructive evaluation of a sample 
using a computer, a display and a computer program imple- 
menting an algorithm, said algorithm includes determining 
the velocity of terahertz electromagnetic radiation in said 
sample without prior knowledge of the thickness of said 
sample and determining said thickness of said sample without 
prior knowledge of said velocity of said terahertz electromag- 
netic radiation in said sample, said terahertz electromagnetic 
radiation produced by a source, namely, a transceiver, spaced 
apart from said sample under test and propagated at the speed 
of light, c, in a medium located between said source, namely, 
said transceiver, and said sample as claimed in claim 1, 
wherein the step of conditioning said front surface (FS) echo 
waveform includes processes selected from the group of 
wavelet denoising, wavelet reconstruction, and filtering. 

6. A process for non-destructive evaluation of a sample 
using a computer, a display and a computer program imple- 
menting an algorithm, said algorithm includes determining 
the velocity of terahertz electromagnetic radiation in said 
sample without prior knowledge of the thickness of said 
sample and determining said thickness of said sample without 
prior knowledge of said velocity of said terahertz electromag- 
netic radiation in said sample, said terahertz electromagnetic 
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radiation produced by a source, namely, a transceiver, spaced 
apart from said sample under test and propagated at the speed 
of light, c, in a medium located between said source, namely, 
said transceiver, and said sample as claimed in claim 1, 
5 wherein the step of determining and storing according to an 
algorithm, on a scan point by scan point basis, the velocity of 
terahertz electromagnetic radiation in said sample without 
prior knowledge of said thickness, is performed according to 
said algorithm solving for V using the equation: V=c(l-At/ 
to 2x) and the step of determining according to said algorithm 
and storing, on a scan point by scan point basis, said thickness 
of said sample without prior knowledge of said velocity of 
said terahertz electromagnetic radiation in said sample, is 
performed according to said algorithm solving for d using the 
15 equation: d=c(2x-At)/2. 

7. A process for non-destructive evaluation of a sample 
using a computer, a display and a computer program imple- 
menting an algorithm, said algorithm includes determining 
the velocity of terahertz electromagnetic radiation in said 

20 sample without prior knowledge of the thickness of said 
sample and determining said thickness of said sample without 
prior knowledge of said velocity of said terahertz electromag- 
netic radiation in said sample, said terahertz electromagnetic 
radiation produced by a source, namely, a transceiver, spaced 
25 apart from said sample under test and propagated at the speed 
of light, c, in a medium located between said source, namely, 
said transceiver, and said sample as claimed in claim 1, 
wherein the step of calculating the difference in time, 2t, 
between said front surface (FS) echo waveforms and said 
30 back surface (BS) echo waveforms on a scan point by scan 
point basis generating and storing a 2t file and the step of 
calculating according to said algorithm the difference in time, 
At, between said front surface (BS) echo waveforms and said 
substrate (M M ) echo waveforms on a scan point by scan point 
35 basis generating and storing a At file is performed according 
to said algorithm by cross-correlating said waveforms on a 
scan point by scan point basis. 

8. A process for non-destructive evaluation of a sample 
using a computer, a display and a computer program imple- 

40 menting an algorithm, said algorithm includes determining 
the velocity of terahertz electromagnetic radiation in said 
sample without prior knowledge of the thickness of said 
sample and determining said thickness of said sample without 
prior knowledge of said velocity of said terahertz electromag- 
45 netic radiation in said sample, said terahertz electromagnetic 
radiation produced by a source, namely, a transceiver, spaced 
apart from said sample under test and propagated at the speed 
of light, c, in a medium located between said source, namely, 
said transceiver, and said sample as claimed in claim 1, 
50 wherein the step of calculating the difference in time, 2 t, 
between said front surface (FS) echo waveforms and said 
back surface (BS) echo waveforms on a scan point by scan 
point basis generating and storing a 2xfile and the step of 
calculating according to said algorithm the difference in time, 
55 At, between said front surface (BS) echo waveforms and said 
substrate (M M ) echo waveforms on a scan point by scan point 
basis generating and storing a At file is performed by com- 
paring the peaks of said waveforms on a scan point by scan 
point basis. 

60 9. A process for non-destructive evaluation of a sample 

using a computer, a display and a computer program imple- 
menting an algorithm, said algorithm includes determining 
the velocity of terahertz electromagnetic radiation in said 
sample without prior knowledge of the thickness of said 
65 sample and determining said thickness of said sample without 
prior knowledge of said velocity of said terahertz electromag- 
netic radiation in said sample, said terahertz electromagnetic 
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radiation produced by a source, namely, a transceiver, spaced 
apart from said sample under test and propagated at the speed 
of light, c, in a medium located between said source, namely, 
said transceiver, and said sample as claimed in claim 1, com- 
prising the steps of: 

wherein said step of fusing on a point by point basis 
includes fusing said gated substrate echo surface (M M ) 
waveform and said gated front surface echo (FS) wave- 
form into said back surface echo (BS) waveform and 
doubling the time base of said fused data waveform. 

10. A process for non-destructive evaluation of a sample 
using a computer, a display and a computer program imple- 
menting an algorithm, said algorithm includes determining 
the velocity of terahertz electromagnetic radiation in said 
sample without prior knowledge of the thickness of said 
sample and determining said thickness of said sample without 
prior knowledge of said velocity of said terahertz electromag- 
netic radiation in said sample, said terahertz electromagnetic 
radiation produced by a source, namely, a transceiver, spaced 
apart from said sample under test and propagated at the speed 
of light, c, in a medium located between said source, namely, 
said transceiver, and said sample, comprising the steps of: 

emitting and scanning terahertz electromagnetic radiation 
from said source, namely, said transceiver, said terahertz 
radiation applied to a plurality of scan points arranged 
according to x and y coordinates of said source, namely, 
said transceiver, said plurality of scan points includes an 
area at least as large as said surface of said sample; 

forming a substrate (M M ) echo waveform file for each scan 
point; 

placing said sample on said substrate; 

forming a front surface (FS) echo waveform file and a back 
surface (BS) echo waveform file for each scan point; 

superimposing, using a graphical user interface, said front 
surface (FS) echo waveform files, saidback surface (BS) 
echo waveform files and said substrate echo (M M ) wave- 
form for each scan point; 

outputting and displaying an image of one of said back 
surface (BS) echo files for inspection of said superim- 
posed files; 

interactively, through a graphical user interface, inspecting 
said superimposed files; and selecting one of said super- 
imposed files for gating; 

gating, using a graphical user interface, said front surface 
(FS) echo waveform and said substrate echo (M M ) wave- 
form and fusing said front surface (FS) echo waveform 
and said substrate echo (M M ) waveform into said back 
surface (BS) echo waveform on a scan point by scan 
point basis and outputting an image of said fused wave- 
forms; 

selecting, using a graphical user interface, one of said fused 
waveforms and gating said waveform for determining 2x 
and At on a scan point by scan point basis; 

determining and storing, based on 2r and At, on a scan 
point by scan point basis, the velocity of terahertz elec- 
tromagnetic radiation in said sample without prior 
knowledge of the thickness of said sample; and, 

determining and storing, based on 2% and At, on a scan 
point by scan point basis, said thickness of said sample 
without prior knowledge of said velocity of said tera- 
hertz electromagnetic radiation in said sample. 

11. A process for non-destructive evaluation of a sample 
using a computer, a display and a computer program imple- 
menting an algorithm, said algorithm includes determining 
the velocity of terahertz electromagnetic radiation in said 
sample without prior knowledge of the thickness of said 
sample and determining said thickness of said sample without 
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prior knowledge of said velocity of said terahertz electromag- 
netic radiation in said sample, said terahertz electromagnetic 
radiation produced by a source, namely, a transceiver, spaced 
apart from said sample under test and propagated at the speed 
5 of light, c, in a medium located between said source, namely, 
said transceiver, and said sample as claimed in claim 10, 
further comprising the steps of: 

displaying, on a scan point by scan point basis, the velocity 
of terahertz electromagnetic radiation in said sample 
to without prior knowledge of the thickness of said sample; 

and, 

displaying, on a scan point by scan point basis, said thick- 
ness of said sample without prior knowledge of said 
velocity of said terahertz electromagnetic radiation in 
15 said sample. 

12. A process for non-destructive evaluation of a sample 
using a computer, a display and a computer program imple- 
menting an algorithm, said algorithm includes determining 
the velocity of terahertz electromagnetic radiation in said 

20 sample without prior knowledge of the thickness of said 
sample and determining said thickness of said sample without 
prior knowledge of said velocity of said terahertz electromag- 
netic radiation in said sample, said terahertz electromagnetic 
radiation produced by a source, namely, a transceiver, spaced 
25 apart from said sample under test and propagated at the speed 
of light, c, in a medium located between said source, namely, 
said transceiver, and said sample as claimed in claim 10, 
further comprising the steps of: 

calculating, on a scan point by scan point basis, the density 
30 of terahertz electromagnetic radiation in said sample 
without prior knowledge of the thickness of said sample. 

13. A process for non- destructive evaluation of a sample 
using a computer, a display and a computer program imple- 
menting an algorithm, said algorithm includes determining 

35 the velocity of terahertz electromagnetic radiation in said 
sample without prior knowledge of the thickness of said 
sample and determining said thickness of said sample without 
prior knowledge of said velocity of said terahertz electromag- 
netic radiation in said sample, said terahertz electromagnetic 
40 radiation produced by a source, namely, a transceiver, spaced 
apart from said sample under test and propagated at the speed 
of light, c, in a medium located between said source, namely, 
said transceiver, and said sample as claimed in claim 12, 
further comprising the steps of: 

45 displaying, on a scan point by scan point basis, the density 
of terahertz electromagnetic radiation in said sample 
without prior knowledge of the thickness of said sample. 

14. A process for non-destructive evaluation of a sample 
using a computer, a display and a computer program imple- 

50 menting an algorithm, said algorithm includes determining 
the velocity of terahertz electromagnetic radiation in said 
sample without prior knowledge of the thickness of said 
sample and determining said thickness of said sample without 
prior knowledge of said velocity of said terahertz electromag- 
55 netic radiation in said sample, said terahertz electromagnetic 
radiation produced by a source, namely, a transceiver, spaced 
apart from said sample under test and propagated at the speed 
of light, c, in a medium located between said source, namely, 
said transceiver, and said sample as claimed in claim 10, 
60 comprising the steps of: 

wherein said step of fusing on a point by point basis 
includes fusing said gated substrate echo surface (M M ) 
waveform and said gated front surface echo (FS) wave- 
form into said back surface echo (BS) waveform and 
65 doubling the time base of said fused data waveform. 

15. A process for non-destructive evaluation of a sample 
using a computer, a computer display and a computer pro- 
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gram implementing an algorithm, said algorithm includes 
determining the velocity of terahertz electromagnetic radia- 
tion in said sample without prior knowledge of the thickness 
of said sample and determining said thickness of said sample 
without prior knowledge of said velocity of said terahertz 
electromagnetic radiation in said sample, said terahertz elec- 
tromagnetic radiation produced by a source, namely, said 
transceiver, spaced apart from said sample under test and 
propagated at the speed of light, c, in a medium located 
between said source, namely, said transceiver, and said 
sample, comprising the steps of: 

emitting and scanning terahertz electromagnetic radiation 
from said source, namely, said transceiver, said terahertz 
radiation applied to a plurality of scan points arranged 
according to x and y coordinates of said source, namely, 
said transceiver, said plurality of scan points includes an 
area at least as large as said surface of said sample; 
recording and storing in a substrate echo waveform data set 
file, on a scan point by scan point basis, said substrate 
echo (M M ) waveform of said terahertz electromagnetic 
radiation received from said substrate without said 
sample present; 

placing said sample on said substrate; 
recording and storing in a back surface echo waveform data 
set file, on a scan point by scan point basis, the back 
surface echo (BS) waveform of said terahertz electro- 
magnetic radiation received from said substrate with 
said sample present; 

recording and storing in a front surface echo waveform 
data set file, on a scan point by scan point basis, the front 
surface echo (FS) waveform of said terahertz electro- 
magnetic radiation received from said front surface (FS) 
of the sample; 

selecting and superimposing, according to said algorithm, 
said front surface (FS) echo waveform of said selected 
scan point, said back surface (BS) echo waveform of 
said selected scan point, and said substrate echo (M M ) 
waveform of said selected scan point; 
gating, interactively, and fusing according to said algo- 
rithm said substrate echo (M") and said front surface 
(FS) echo waveform forming a fused waveform on a 
scan point by scan point basis; 
selecting, interactively, one of said fused waveforms and 
gating said waveform for determining 2x and At on a 
scan point by scan point basis; 
determining, storing, and displaying, based on 2x and At, 
on a scan point by scan point basis, the velocity of 
terahertz electromagnetic radiation in said sample with- 
out prior knowledge of the thickness of said sample, said 
velocity, V, determined by solving the equation, V=c(l - 
At/2x); and, 

determining, storing, and displaying, based on 2x and At, 
on a scan point by scan point basis, said thickness of said 
sample without prior knowledge of said velocity of said 
terahertz electromagnetic radiation in said sample, said 
thickness, d, determined by solving the equation, d=c 
(2x-At)/2. 

16. A process for non-destructive evaluation of a sample 
using a computer, a computer display and a computer pro- 
gram implementing an algorithm, said algorithm includes 
determining the velocity of terahertz electromagnetic radia- 
tion in said sample without prior knowledge of the thickness 
of said sample and determining said thickness of said sample 
without prior knowledge of said velocity of said terahertz 
electromagnetic radiation in said sample, said terahertz elec- 
tromagnetic radiation produced by a source spaced apart from 
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said sample under test and propagated at the speed of light, c, 
in a medium located between said source and said sample, 
comprising the steps of: 

emitting and scanning terahertz electromagnetic radiation 
5 from said source said terahertz radiation applied to a 
plurality of scan points arranged according to x and y 
coordinates of said scanning, said plurality of scan 
points includes an area at least as large as said surface of 
said sample; 

l o recording and storing in a sub strate echo waveform data set 

file, on a scan point by scan point basis, a substrate echo 
(M") waveform of said terahertz electromagnetic radia- 
tion received from said substrate without said sample 
present, and recording and storing a front surface (FS) 
15 waveform and back surface (BS) waveform of said tera- 

hertz electromagnetic radiation received with the 
sample present; 

superimposing said substrate echo (M M ) waveform, said 
front surface (FS) echo waveform, and said back surface 
20 (BS) echo waveform together; 

gating and conditioning said substrate echo (M") and said 
front surface (FS) echo waveform; 

fusing said substrate echo (M M ) waveform and said front 
surface echo (FS) waveform into said back surface (BS) 
25 echo waveform; 

gating said echo waveforms and determining said time 
difference between said front surface echo (FS) wave- 
form and said back surface (BS) echo waveform, 2x, and 
storing values of 2x for each scan point; 

30 gating said echo waveforms and determining said time 
difference between said substrate echo wave (M M ) wave- 
form and said back surface (BS) waveform, At, and 
storing values of At for each scan point; 

using said stored values of 2x and At in said algorithm to 
35 determine the velocity of the terahertz electromagnetic 

radiation in said sample on a scan point by scan point 
basis, said algorithm includes the step of determining 
and storing, on a scan point by scan point basis, the 
velocity of terahertz electromagnetic radiation in said 
40 sample without prior knowledge of said thickness, said 
velocity, V, is determined by solving the equation, V=c 
(l-At/2x), for each said scan point and storing said 
velocity value in a velocity computer file on a scan point 
by scan point basis; 

45 using said stored values of 2x and At in said algorithm to 
determine the thickness of said sample on a scan point by 
scan point basis, said algorithm includes the step of 
determining on a scan point by scan point basis, said 
thickness, d, of said sample without prior knowledge of 
50 said velocity of said terahertz electromagnetic radiation 
in said sample, is determined by solving the equation, 
d=c(2x-At)/2, for each said scan point and storing said 
thickness value in a thickness computer file. 

17. A process for non-destructive evaluation of a sample 
55 using a computer, a computer display and a computer pro- 
gram implementing an algorithm, said algorithm includes 
determining the velocity of terahertz electromagnetic radia- 
tion in said sample without prior knowledge of the thickness 
of said sample and determining said thickness of said sample 
60 without prior knowledge of said velocity of said terahertz 
electromagnetic radiation in said sample, said terahertz elec- 
tromagnetic radiation produced by a source spaced apart from 
said sample under test and propagated at the speed of light, c, 
in a medium located between said source and said sample, as 
65 claimed in claim 16, wherein the steps of gating and condi- 
tioning said substrate echo (M") and said front surface (FS) 
echo waveform include delaying said substrate echo in time 



US 7,933,027 B1 


47 

to better gate said substrate echo and then remove said delay 
when said step of fusing said substrate echo (M M ) waveform 
and said front side echo (FS) waveform into said back side 
echo waveform (BS) is performed. 

18. A process for non-destructive evaluation of a sample 
using a computer, a computer display and a computer pro- 
gram implementing an algorithm, said algorithm includes 
determining the velocity of terahertz electromagnetic radia- 
tion in said sample without prior knowledge of the thickness 
of said sample and determining said thickness of said sample 
without prior knowledge of said velocity of said terahertz 
electromagnetic radiation in said sample, said terahertz elec- 
tromagnetic radiation produced by a source spaced apart from 
said sample under test and propagated at the speed of light, c, 
in a medium located between said source and said sample, as 
claimed in claim 17, wherein said steps of: gating said echo 
waveforms and determining said time difference between 
said front surface echo (FS) waveform and said back surface 
(BS) echo waveform, 2x, and storing values of 2x for each 
scan point; and, gating said echo waveforms and determining 
said time difference between said substrate echo wave (M M ) 
waveform and said back surface (BS) waveform, At, and 
storing values of At are performed by said algorithm cross- 
correlating said waveforms for each scan point. 
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19. A process for non-destructive evaluation of a sample 
using a computer, a computer display and a computer pro- 
gram implementing an algorithm, said algorithm includes 
determining the velocity of terahertz electromagnetic radia- 
5 tion in said sample without prior knowledge of the thickness 
of said sample and determining said thickness of said sample 
without prior knowledge of said velocity of said terahertz 
electromagnetic radiation in said sample, said terahertz elec- 
tromagnetic radiation produced by a source spaced apart from 
l o said sample under test and propagated at the speed of light, c , 
in a medium located between said source and said sample, as 
claimed in claim 17, wherein said steps of: gating said echo 
waveforms and determining said time difference between 
said front surface echo (FS) waveform and said back surface 
15 (BS) echo waveform, 2x, and storing values of 2x for each 
scan point; and, gating said echo waveforms and determining 
said time difference between said substrate echo wave (M M ) 
waveform and said back surface (BS) waveform, At, and 
storing values of At are performed by said algorithm by com- 
20 parison of peak to peak amplitudes for each scan point. 



